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We use the macroscopic quantum electrodynamics approach suitable for ab- 
sorbing and dispersing media to study the properties and role of collective 
surface excitations — excitons and plasmons — in single-wall and double- 
wall carbon nanotubes. We show that the interactions of excitonic states 
with surface electromagnetic modes in individual small-diameter ( < 1 nm) 
single-walled carbon nanotubes can result in strong exciton-surface-plasmon 
coupling. Optical response of individual nanotubes exhibits Rabi splitting 
~0.1 eV, both in the linear excitation regime and in the non-linear excitation 
regime with the photoinduced biexcitonic states formation, as the exciton 
energy is tuned to the nearest interband surface plasmon resonance of the 
nanotube. An electrostatic field applied perpendicular to the nanotube axis 
can be used to control the exciton-plasmon coupling. For double-wall carbon 
nanotubes, we show that at tube separations similar to their equilibrium dis- 
tances interband surface plasmons have a profound effect on the inter-tube 
Casimir force. Strong overlapping plasmon resonances from both tubes war- 
rant their stronger attraction. Nanotube chiralities possessing such collective 
excitation features will result in forming the most favorable inner-outer tube 
combination in double- wall carbon nanotubes. These results pave the way 
for the development of new generation of tunable optoelectronic and nano- 
electromechanical device applications with carbon nanotubes. 
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1 Introduction 

Single- walled carbon nanotubes (CNs) are quasi-one-dimensional (ID) cylindrical 
wires consisting of graphene sheets rolled-up into cylinders with diameters ~ 1 — 
10 nm and lengths ~ 1 — 10^ [H [2l [3l S]. CNs are shown to be useful as 
miniaturized electronic, electromechanical, and chemical devices [5], scanning probe 
devices [6], and nanomaterials for macroscopic composites [7]. The area of their 
potential applications was recently expanded to nanophotonics [8l |9l |10l [TTl [121 E] 
after the demonstration of controllable single-atom incapsulation into CNs |14^ [T5| 
[T6| [T7] . and even to quantum cryptography since the experimental evidence was 
reported for quantum correlations in the photoluminescence spectra of individual 
nanotubes fl8\. 

For pristine (undoped) single-walled CNs, the numerical calculations predicting 
large exciton binding energies (~ 0.3—0.6 eV) in semiconducting CNs |19 1 [20 l [2T] and 
even in some small-diameter (~0.5 nm) metallic CNs [22], followed by the results 
of various exciton photoluminescence measurements [El [231 EH [251 EHl [27] , have 
become available. These works, together with other reports investigating the role 
of effects such as intrinsic defects [25l [28] , exciton-phonon interactions [261 ESI [291 
[30l l3Tj , biexciton formation [32l [33] , exciton-surface-plasmon coupling [SH [351 IMl 
[37] . external magnetic [381 El] and electric fields [371 HO], reveal the variety and 
complexity of the intrinsic optical properties of CNs [41] . 

Carbon nanotubes combine advantages such as electrical conductivity, chemical 
stability, high surface area, and unique optoelectronic properties that make them 
excellent potential candidates for a variety of applications, including efficient solar 
energy conversion [7], energy storage [14J, optical nanobiosensorics [32]. However, 
the quantum yield of individual CNs is normally very low. Nanotube composites 
of CN bundles and/or films could surpass this difficulty, opening up new paths for 
the development of high-yield, high-performance optoelectronics applications with 
CNs |43[ 144). Understanding the inter-tube interactions is important in order to be 
able to tailor properties of CN bundles and films, as well as properties of multi-wall 
CNs. This is also important for experimental realization of new effects and devices 
proposed recently, such as trapping of cold atoms [42[I45] and their entanglement [TT] 
near single- walled CNs, surface profiling [6] and nanolithography applications [l6] 
with double-wall CNs. 

Here, we use the macroscopic Quantum ElectroDynamics (QED) formalism de- 
veloped earlier for absorbing and dispersive media [371 HSl [Ml [9] and then success- 
fully employed to study near-field EM effects in hybrid CN systems [lOl [TTl US], 
to investigate the properties and role of collective surface excitations — excitons 
and plasmons — in single- wall and double- wall CNs. First, we show that, due 
to the presence of low-energy (~ 0.5 — 2 eV) weakly-dispersive interband plasmon 
modes [50] and large exciton excitation energies in the same energy domain |51[ [52] , 
the excitons can form strongly coupled mixed exciton-plasmon excitations in indi- 
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vidual small-diameter (< 1 nm) semiconducting single-walled CNs. The exciton- 
plasmon coupling (and the exciton emission accordingly) can be controlled by an 
external electrostatic field applied perpendicular to the CN axis (the quantum con- 
fined Stark effect) . The optical response of individual CNs exhibits the Rabi splitting 
of ~ 0.1 eV, both in the linear excitation regime and in the non- linear excitation 
regime with the photoinduced biexcitonic states formation, as the exciton energy is 
tuned to the nearest interband surface plasmon resonance of the CN. Previous stud- 
ies of the exciton-plasmon coupling have been focused on artificially fabricated hy- 
brid plasmonic nanostructures, such as dye molecules in organic polymers deposited 
on metallic films [53], semiconductor quantum dots coupled to metallic nanoparti- 
cles [54], or nanowires [55]. where semiconductor material carries the exciton and 
metal carries the plasmon. Our results are particularly interesting since they reveal 
the fundamental electromagnetic (EM) phenomenon — the strong exciton-plasmon 
coupling — in an individual quasi- ID nanostructure, a carbon nanotube, as well as 
its tunability feature by means of the quantum confined Stark effect. We expect 
these results to open up new paths for the development of tunable optoelectronic 
device applications with optically excited carbon nanotubes, including the strong 
excitation regime with optical non-linearities. 

Next, we turn to the double- wall carbon nanotubes to investigate the effect of 
collective surface excitations on the inter-tube Casimir interaction in these systems. 
The Casimir interaction is a paradigm for a force induced by quantum EM fiuctua- 
tions. The fundamental nature of this force has been studied extensively ever since 
the prediction of the existence of an attraction between neutral metallic mirrors in 
vacuum |491 156] . In recent years, the Casimir effect has acquired a much broader 
impact due to its importance for nanostructured materials and devices. The de- 
velopment and operation of micro- and nano-electromechanical systems are limited 
due to unwanted effects, such as stiction, friction, and adhesion, originating from the 
Casimir force |57] . This interaction is also an important component for the stability 
of nanomaterials. Here, we show that at tube separations similar to their equilib- 
rium distances interband surface plasmons have a profound effect on the inter-tube 
Casimir force. Strong overlapping plasmon resonances from both tubes warrant 
their stronger attraction. Nanotube chiralities possessing such collective excitation 
features will result in forming the most favorable inner-outer tube combination in 
double-wall carbon nanotubes. This theoretical understanding is important for the 
development of nano-electromechanical devices with CNs. 

This Chapter is organized as follows. Section [2] introduces the general Hamilto- 
nian of the exciton interaction with vacuum-type quantized surface EM modes of 
a single-walled CN. No external EM field is assumed to be applied. The vacuum- 
type-field we consider is created by CN surface EM fluctuations. Section [3] explains 
how the interaction introduced results in the coupling of the excitonic states to 
the nanotube's surface plasmon modes. Here we derive and discuss the character- 
istics of the coupled exciton-plasmon excitations, such as the dispersion relation. 
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Figure 1: The geometry of the problem. 



the plasmon density of states (DOS), and the optical response functions, for par- 
ticular semiconducting CNs of different diameters. We also analyze how the elec- 
trostatic field applied perpendicular to the CN axis affects the CN band gap, the 
exciton binding energy, and the surface plasmon energy, to explore the tunability 
of the exciton-surface-plasmon coupling in CNs. Section [4] derives and analyzes the 
Casimir interaction between two concentric cylindrical graphene sheets comprising 
a double-wall CN. The summary and conclusions of the work are given in Sec. [5l All 
the technical details about the construction and diagonalization of the exciton-field 
Hamiltonian, the EM field Green tensor derivation, the perpendicular electrostatic 
field effect, are presented in the Appendices in order not to interrupt the flow of the 
arguments and results. 

2 Exciton-electromagnetic-field interaction 
on the nanotube surface 

We consider the vacuum-type EM interaction of an exciton with the quantized 
surface electromagnetic fluctuations of a single-walled semiconducting CN by using 
our recently developed Green function formalism to quantize the EM field in the 
presence of quasi- ID absorbing bodies [581 EHl EOl Ell ESI [9]. No external EM field 
is assumed to be applied. The nanotube is modelled by an infinitely thin, infinitely 
long, anisotropically conducting cylinder with its surface conductivity obtained from 
the realistic band structure of a particular CN. Since the problem has the cylindrical 
symmetry, the orthonormal cylindrical basis {er,e^,ez} is used with the vector 
directed along the nanotube axis as shown in Fig.[Tl Only the axial conductivity, azz, 
is taken into account, whereas the azimuthal one, a^y,, being strongly suppressed 
by the transverse depolarization effect [Ml EH ESI ESI ES [68] , is neglected. 
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The total Hamiltonian of the coupled exciton-photon system on the nanotube 
surface is of the form 

H = Hp + Hex + Hint, (1) 

where the three terms represent the free (medium-assisted) EM field, the free (non- 
interacting) exciton, and their interaction, respectively. More explicitly, the second 
quantized field Hamiltonian is 

HF = y^ diohup{n,L,)f{n,u;), (2) 

where the scalar bosonic field operators p{n,uj) and /(n, cj) create and annihilate, 
respectively, the surface EM excitation of frequency u at an arbitrary point n = 
Rn = {RcN,'^n,Zn} associated with a carbon atom (representing a lattice site - 
Fig. [1]) on the surface of the CN of radius RcN- The summation is made over all the 
carbon atoms, and in the following it is replaced by the integration over the entire 
nanotube surface according to the rule 

^ dR„... = — / dipnRcN dzn..., (3) 
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where Sq = {3V3/4)b'^ is the area of an elementary equilateral triangle selected 
around each carbon atom in a way to cover the entire surface of the nanotube, 
6=1.42 A is the carbon-carbon interatomic distance. 

The second quantized Hamiltonian of the free exciton (see, e.g., Ref. [69]) on the 
CN surface is of the form 

Hex = J2Ef{n)Bl^^jB^j = ^Ef{k)BlfB^j, (4) 

n,mj k,/ 

where the operators B^ j- and Bnj create and annihilate, respectively, an exciton 
with the energy Ef{n) in the lattice site n of the CN surface. The index / (7^ 0) 
refers to the internal degrees of freedom of the exciton. Alternatively, 

= E</«^'" B^j = {Blj)^ (5) 

create and annihilate the /-internal-state exciton with the quasi-momentum k = 
{k^,kz}, where the azimuthal component is quantized due to the transverse con- 
finement effect and the longitudinal one is continuous, is the total number of the 
lattice sites (carbon atoms) on the CN surface. The exciton total energy is then 
written in the form 

k2h.2 
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Here, the first term represents the excitation energy 

Eglik^) = Eg{k^) + El^\k^) (7) 

of the /-internal-state exciton with the (negative) binding energy E^ , created via 
the interband transition with the band gap 

Eg{k^) = £e{k^) + £hik^), (8) 

where £e^h are transversely quantized azimuthal electron-hole subbands (see the 
schematic in Fig. [2|). The second term in Eq. ([6|) represents the kinetic energy of the 
translational longitudinal movement of the exciton with the effective mass M^x = 
iTie + mfi, where me and ruh are the (subband-dependent) electron and hole effective 
masses, respectively. The two equivalent free-exciton Hamiltonian representations 
are related to one another via the obvious orthogonality relationships 

n k 

with the k-summation running over the first Brillouin zone of the nanotube. The 
bosonic field operators in Hp are transformed to the k-representation in the same 
way. 

The most general (non-relativistic, electric dipole) exciton-photon interaction on 
the nanotube surface can be written in the form (we use the Gaussian system of 
units and the Coulomb gauge; see details in Appendix A) 

duj [g'j^\n,m,uj)Blj - gj"^(n, m, a;)5nj ] /(m, + /i.e., (10) 

n,m,/ - 

where 



'•'■^^''n, m, w) = gi (n, m, w) ± — g'|(n, m, lo) (11) 



with 



,/'l)(n,m,c^) = -i^^T,hw^^a,,{RcN.^) (d^). ^("^^..(n, m, w) (12) 

being the interaction matrix element where the exciton with the energy -Eexc = ^/ 
is excited through the electric dipole transition (d;^)^ = (OKdn)^!/) in the lattice site 
n by the nanotube's transversely (longitudinally) polarized surface EM modes. The 
modes are represented in the matrix element by the transverse (longitudinal) part of 
the Green tensor zz-component Gj;2(n, m, w) of the EM subsystem (Appendix B). 
This is the only Green tensor component we have to take into account. All the other 
components can be safely neglected as they are greatly suppressed by the strong 
transverse depolarization effect in CNs [63l [Ml ESI IMl EZl [68] • As a consequence, 
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Figure 2: Schematic of the two transversely quantized azimuthal electron-hole sub- 
bands {left), and the first-interband ground-internal-state exciton energy [right) in 
a small-diameter semiconducting carbon nanotube. Subbands with indices j = 1 
and 2 are shown, along with the optically allowed (exciton-related) interband tran- 
sitions [67J. See text for notations. 

only azz{RcN,^), the axial dynamic surface conductivity per unit length, is present 
in Eq.(II2]). 

Equations p!])- (|12p form the complete set of equations describing the exciton- 
photon coupled system on the CN surface in terms of the EM field Green tensor 
and the CN surface axial conductivity. 

3 Exciton-surface-plasmon coupling 

For the following it is important to realize that the transversely polarized sur- 
face EM mode contribution to the interaction ()10p - ()12p is negligible compared to 
the longitudinally polarized surface EM mode contribution. As a matter of fact, 
"'"GzzCn, m, w) = in the model of an infinitely thin cylinder we use here (Ap- 
pendix B), thus yielding 

g|(n,m,a;)=0, g'^\n,in,uj) = ±— g\{n,m, uj) (13) 

in Eqs. p^ ~ p^ . The point is that, because of the nanotube quasi-one-dimen- 
sionality, the exciton quasi-momentum vector and all the relevant vectorial matrix 
elements of the momentum and dipole moment operators are directed predominantly 
along the CN axis (the longitudinal exciton; see, however, Ref. [70]). This prevents 
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the exciton from the electric dipole couphng to transversely polarized surface EM 
modes as they propagate predominantly along the CN axis with their electric vectors 
orthogonal to the propagation direction. The longitudinally polarized surface EM 
modes are generated by the electronic Coulomb potential (see, e.g., Ref. [71]), and 
therefore represent the CN surface plasmon excitations. These have their electric 
vectors directed along the propagation direction. They do couple to the longitudi- 
nal excitons on the CN surface. Such modes were observed in Ref. [50j . They occur 
in CNs both at high energies (well-known 7r-plasmon at ~ 6 eV) and at compara- 
tively low energies of ~ 0.5 — 2 eV. The latter ones are related to the transversely 
quantized interband (inter-van Hove) electronic transitions. These weakly-dispersive 
modes [501172] are similar to the intersubband plasmons in quantum wells [73]. They 
occur in the same energy range of ~ 1 eV where the exciton excitation energies are 
located in small-diameter (< 1 nm) semiconducting CNs [5H I52j . In what follows 
we focus our consideration on the exciton interactions with these particular surface 
plasmon modes. 



3.1 The dispersion relation 

To obtain the dispersion relation of the coupled exciton-plasmon excitations, we 
transfer the total Hamiltonian (fT])- pop and (I13p to the k-representation using 
Eqs. ([5]) and Q, and then diagonalize it exactly by means of Bogoliubov's canonical 
transformation technique (see, e.g., Ref. [74]). The details of the procedure are given 
in Appendix C. The Hamiltonian takes the form 

H=J2 ^M(k) el(k)4(k) + Eo . (14) 

k,M=l,2 

Here, the new operator 

4(k) = [K{K^f)B^j - v,{k,uj)Bi^j] (15) 
/ 



+ jlu; [u^ik,u)f{k,u;)-v;{k,u;)f\-k,u;) 



annihilates and ^jj(k) = [^^(k)]^^ creates the exciton-plasmon excitation of branch 
/i, the quantities and are appropriately chosen canonical transformation co- 
efficients. The "vacuum" energy Eq represents the state with no exciton-plasmons 
excited in the system, and ^^(k) is the exciton-plasmon energy given by the solu- 
tion of the following (dimensionless) dispersion relation 

2 2 n xTl(x)p(x) , , 

dx—^^ = 0. (16) 

^ Jo Xfj. ~ X 

Here, 

x = ^, x, = ^, e, = ^ (17) 
270 270 270 
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with 7o = 2.7 eV being the carbon nearest neighbor overlap integral entering the CN 
surface axial conductivity azz{RcN ,^)- The function 



with di = Yl,n{^\{^n)z\f) represents the (dimensionless) spontaneous decay rate, and 

p[x) = ^-^^ Re^-- (19) 

stands for the surface plasmon density of states (DOS) which is responsible for the 
exciton decay rate variation due to its coupling to the plasmon modes. Here, a = 
e^/^c= 1/137 is the fine-structure constant and (Tzz = '^'^^(^zz/ is the dimensionless 
CN surface axial conductivity per unit length. 

Note that the conductivity factor in Eq. (fT9]) equals 



1 4ac / n \ 1 

Re _ , , = Im — (20) 

(yzz(.x) RcN y^ioxj ezz{x) - 1 



in view of Eq. (fT7|) and equation 



lUJ 



(Tzz{x) = --^[ezz{x)-l\ (21) 

representing the Drude relation for CNs, where is the longitudinal (along the 
CN axis) dielectric function, S and pT are the surface area of the tubule and the 
number of tubules per unit volume, respectively |59 1I621 [M]. This relates very closely 
the surface plasmon DOS function (jlOp to the loss function — Im(l/e) measured in 
Electron Energy Loss Spectroscopy (EELS) experiments to determine the properties 
of collective electronic excitations in solids |50j . 

Figure [3] shows the low-energy behaviors of azz{x) and Re[l/c722(x)] for the 
(11,0) and (10,0) CNs {RcN = 0.43 nm and 0.39 nm, respectively) we study here. 
We obtained them numerically as follows. First, we adapt the nearest-neighbor non- 
orthogonal tight-binding approach [75j to determine the realistic band structure of 
each CN. Then, the room-temperature longitudinal dielectric functions Czz are cal- 
culated within the random-phase approximation j76t [77| , which are then converted 
into the conductivities (Tzz by means of the Drude relation. Electronic dissipation 
processes are included in our calculations within the relaxation-time approximation 
(electron scattering length of 12,QRcn was used [30]). We did not include excitonic 
many-electron correlations, however, as they mostly affect the real conductivity 
Re(a^^) which is responsible for the CN optical absorption \22 \ Wl\. whereas we 
are interested here in Re(l/a^2) representing the surface plasmon DOS according to 
Eq. (jl9p . This function is only non-zero when the two conditions, Im[a'22(a;)] = 
and Re[fT22(x)] — )• 0, are fulfilled simultaneously [72l [73l [76]. These result in the 
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Figure 3: (a),(b) Calculated dimensionless (see text) axial surface conductivities for 
the (11,0) and (10,0) CNs. The dimensionless energy is defined as [Energy]/2^Q, 
according to Eq. (fTTl) . 

peak structure of the function Jie{l/azz) as is seen in Fig. [3l It is also seen from 
the comparison of Fig. [3] (b) with Fig. [3] (a) that the peaks broaden as the CN 
diameter decreases. This is consistent with the stronger hybridization effects in 
smaller-diameter CNs [78j. 

Left panels in Figs. H^a) and \Mjo) show the lowest-energy plasmon DOS res- 
onances calculated for the (11,0) and (10,0) CNs as given by the function p{x) 
in Eq. (jlOp . Also shown there are the corresponding fragments of the functions 
^^[^■^^(x)] and Im[a'2^(x)]. In all graphs the lower dimensionless energy limits are 
set up to be equal to the lowest bright exciton excitation energy [E'ixc^ = 1-21 eV 
(x = 0.224) and 1.00 eV {x = 0.185) for the (11,0) and (10,0) CN, respectively, 
as reported in Ref. |51j by directly solving the Bethe-Salpeter equation]. Peaks in 
p{x) are seen to coincide in energy with zeros of Im[CT22(x)] {or zeros of Re[e^^(x)]}, 
clearly indicating the plasmonic nature of the CN surface excitations under con- 
sideration \72\ I79j. They describe the surface plasmon modes associated with the 
transversely quantized interband electronic transitions in CNs [72]. As is seen in 
Fig. m (and in Fig. [3|), the interband plasmon excitations occur in CNs slightly above 
the first bright exciton excitation energy [67], in the frequency domain where the 
imaginary conductivity (or the real dielectric function) changes its sign. This is a 
unique feature of the complex dielectric response function, the consequence of the 
general Kramers-Kronig relation |47j . 

We further take advantage of the sharp peak structure of p(x) and solve the 
dispersion equation ()16p for x^ analytically using the Lorentzian approximation 
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Figure 4: (a),(b) Surface plasmon DOS and conductivities (left panels), and lowest 
bright exciton dispersion when coupled to plasmons (right panels) in (11,0) and 
(10,0) CNs, respectively. The dimensionless energy is defined as [Energy]/2^Q, ac- 
cording to Eq. (jl7p . See text for the dimensionless quasi- momentum. 



Here, Xp and Axp are, respectively, the position and the half-width-at-half-maximum 
of the plasmon resonance closest to the lowest bright exciton excitation energy in 
the same nanotube (as shown in the left panels of Fig. [4]). The integral in Eq. (|16p 
then simplifies to the form 



vr 



xTl{x)p{x) F{xp)Axl 



dx 



F( 



dx 



{x 



arctan 



X. 

+ 2 



|2 + Axj 



with F{xp) = 2xpTQ{xp)p{xp) /t: . This expression is valid for all x^ apart from those 
located in the narrow interval {xp — /S.Xp, Xp + Axp) in the vicinity of the plasmon res- 
onance, provided that the resonance is sharp enough. Then, the dispersion equation 
becomes the biquadratic equation for x^ with the following two positive solutions 
(the dispersion curves) of interest to us 
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£/ + ^ 1 




(23) 



2 2 



Here, = AF(xp)Axp{TT — Axp/xp) with the arctan-function expanded to hnear 
terms in Axp/xp <C 1. 

The dispersion curves (I23p are shown in the right panels in Figs. [Hj^a) and|4ljb) as 
functions of the dimensionless longitudinal quasi-momentum. In these calculations, 
we estimated the interband transition matrix element in r^(xp) [Eq.([l8])] from the 
equation = 3/iA^/4r™'^ according to Hanamura's general theory of the exciton 
radiative decay in spatially confined systems jSOj, where r™"^ is the exciton intrinsic 
radiative lifetime, and A = 2Trch/E with E being the exciton total energy given in 
our case by Eq. ([6]). For zigzag-type CNs considered here, the first Brillouin zone of 
the longitudinal quasi-momentum is given by —2Trh/3b < hkz < 2TTh/3b [T, "2]. The 
total energy of the ground-internal-state exciton can then be written as E = E^xc + 
(27r^/36)^t^/2Mex with — 1 < t < 1 representing the dimensionless longitudinal 
quasi-momentum. In our calculations we used the lowest bright exciton parameters 
eH^J = 1.21 eV and 1.00 eV, r^f = 14.3 ps and 19.1 ps, Mex = 0.44mo and 0.19mo 
(mo is the free-electron mass) for the (11,0) CN and (10,0) CN, respectively, as 
reported in Ref.[5T| by directly solving the Bethe-Salpeter equation. 

Both graphs in the right panels in Fig. |4] are seen to demonstrate a clear an- 
ticrossing behavior with the (Rabi) energy splitting ~ 0.1 eV. This indicates the 
formation of the strongly coupled surface plasmon-exciton excitations in the nan- 
otubes under consideration. It is important to realize that here we deal with the 
strong exciton-plasmon interaction supported by an individual quasi-lD nanostruc- 
ture — a single-walled (small-diameter) semiconducting carbon nanotube, as op- 
posed to the artificially fabricated metal-semiconductor nanostructures studied pre- 
viuosly |53l [54l [55] where the metallic component normally carries the plasmon and 
the semiconducting one carries the exciton. It is also important that the effect comes 
not only from the height but also from the width of the plasmon resonance as it 
is seen from the definition of the Fp factor in Eq. (|23p . In other words, as long as 
the plasmon resonance is sharp enough (which is always the case for interband plas- 
mons), so that the Lorentzian approximation (I22p applies, the effect is determined 
by the area under the plasmon peak in the DOS function (jl9p rather than by the 
peak height as one would expect. 

However, the formation of the strongly coupled exciton-plasmon states is only 
possible if the exciton total energy is in resonance with the energy of a surface plas- 
mon mode. The exciton energy can be tuned to the nearest plasmon resonance in 
ways used for excitons in semiconductor quantum microcavities — thermally \81\ 
[82} [83] (by elevating sample temperature), or /and electrostatically (SH [85l [861 EZ] 
(via the quantum confined Stark effect with an external electrostatic field applied 
perpendicular to the CN axis). As is seen from Eqs. ^ and ([7]), the two possibilities 
influence the different degrees of freedom of the quasi-lD exciton — the (longitudi- 
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nal) kinetic energy and the excitation energy, respectively. Below we study the (less 
trivial) electrostatic field effect on the exciton excitation energy in CNs. 

3.2 The perpendicular electrostatic field effect 

The optical properties of semiconducting CNs in an external electrostatic field di- 
rected along the nanotube axis were studied theoretically in Ref. [40j. Strong oscil- 
lations in the band-to-band absorption and the quadratic Stark shift of the exciton 
absorption peaks with the field increase, as well as the strong field dependence of the 
exciton ionization rate, were predicted for CNs of different diameters and chiralities. 
Here, we focus on the perpendicular electrostatic field orientation. We study how 
the electrostatic field applied perpendicular to the CN axis affects the CN band gap, 
the exciton binding/excitation energy, and the interband surface plasmon energy, 
to explore the tunability of the strong exciton-plasmon coupling effect predicted 
above. The problem is similar to the well-known quantum confined Stark effect first 
observed for the excitons in semiconductor quantum wells [84^ I85j . However, the 
cylindrical surface symmetry of the excitonic states brings new peculiarities to the 
quantum confined Stark effect in CNs. In what follows we will generally be inter- 
ested only in the lowest internal energy (ground) excitonic state, and so the internal 
state index / in Eqs. Q and ([7]) will be omitted for brevity. 

Because the nanotube is modelled by a continuous, infinitely thin, anisotropically 
conducting cylinder in our macroscopic QED approach, the actual local symmetry 
of the excitonic wave function resulted from the graphene Brillouin zone structure is 
disregarded in our model (see, e.g., reviews |4H I67j). The local symmetry is implic- 
itly present in the surface axial conductivity though, which we calculate beforehand 
as described aboveH 

We start with the Schrodinger equation for the electron and hole on the CN 
surface, located at = {RcN,^e, Ze} and = {RcN,^h, Zh}, respectively. They 
interact with each other through the Coulomb potential y(re,r/j) = — e^/e|re — 
r/i|, where e = €zz{0)- The external electrostatic field F = {F, 0,0} is directed 
perpendicular to the CN axis (along the x-axis in Fig. [1]). The Schrodinger equation 
is of the form 



^In real CNs, the existence of two equivalent energy valleys in the 1st Brillouin zone, the K- 
and i^'-valleys with opposite electron helicities about the CN axis, results into dark and bright 
excitonic states in the lowest energy spin-singlet manifold ^8]. Since the electric interaction does 
not involve spin variables, both K- and /(''-valleys are affected equally by the electrostatic field in 
our case, and the detailed structure of the exciton wave function multiplet is not important. This 
is opposite to the non-zero magnetostatic field case where the field affects the K- and if'-valleys 
differently either to brighten the dark excitonic states [35], or to create Landau sublevels [57] for 
longitudinal and perpendicular orientation, respectively. 




(24) 
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with 



2m. 



e,h 



P2 



+ 



9^ 



(25) 



We further separate out the translational and relative degrees of freedom of the 
electron-hole pair by transforming the longitudinal (along the CN axis) motion of 
the pair into its center-of-mass coordinates given hy Z = {rueZe + ruhZh) / Mex and 
z = Ze — Zh- The exciton wave function is approximated as follows 



^'(re,r,,) = e'^^^ (l)ex{z)i^e{^e)i^h{^h)- 



(26) 



The complex exponential describes the exciton center-of-mass motion with the lon- 
gitudinal quasi-momentum along the CN axis. The function (j)ex{z) represents 
the longitudinal relative motion of the electron and the hole inside the exciton. The 
functions tpe{^e) and iphi'^h) are the electron and hole subband wave functions, 
respectively, which represent their confined motion along the circumference of the 
cylindrical nanotube surface. 

Each of the functions is assumed to be normalized to unity. Equations (j24p and 
(j25p are then rewritten in view of Eqs. ©-([HI) to yield 



^2—^-^ + eRcNFcos{iph) 



iph{^h)=£:htph{v>h), 



' 2/i 



+ ^eff(^) 



(t)ex{z) = Eh(f)exiz), 



(27) 

(28) 
(29) 



where fi = memfi/Mex is the exciton reduced mass, and VqQ is the effective longitu- 
dinal electron-hole Coulomb interaction potential given by 



2lT 



2tt 



(30) 



with V being the original electron-hole Coulomb potential written in the cylindrical 
coordinates as 



1 



{z2+4i?2,^sin2[((^,-(^0/2]}i/2' 



(31) 



The exciton problem is now reduced to the ID equation (j29p . where the exciton 
binding energy does depend on the perpendicular electrostatic field through the 
electron and hole subband functions ^jJe,h given by the solutions of Eqs. ([27|) and 
()28p and entering the effective electron-hole Coulomb interaction potential pop . 
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The set of Eqs. (j27p - (|3ip is analyzed in Appendix D. One of the main results 
obtained in there is that the effective Coulomb potential (|30p can be approximated 
by an attractive cusp-type cutoff potential of the form 



e2 



e[\z\+ zo{j,F)] 

where the cutoff parameter zq depends on the perpendicular electrostatic field 
strength and on the electron- hole azimuthal transverse quantization index j = 1,2,... 
(excitons are created in interband transitions involving valence and conduction sub- 
bands with the same quantization index [67] as shown in Fig. [2]) . Specifically, 

7r-21n2[l - A,'(F)1 , , 

with Aj(F) given to the second order approximation in the electric field by 

A,(F) « 2^Me. , (34) 

.... _ ^(i-2) + 1 



l-2j l + 2j' 

where 9{x) is the unit step function. Approximation (j32p is formally valid when 
zo{j, F) is much less than the exciton Bohr radius (= e/i^//ie^) which is estimated 
to be ~ WRcN for the first (j = 1 in our notations here) exciton in CNs [E]. As is 
seen from Eqs. (j33p and (j34p . this is always the case for the first exciton for those 
fields where the perturbation theory applies, i. e. when Ai{F) < 1 in Eq. (I34p . 

Equation ([29]) with the potential (f32|) formally coincides with the one studied by 
Ogawa and Takagahara in their treatments of excitonic effects in ID semiconductors 
with no external electrostatic field applied [5S]. The only diff'erence in our case is 
that our cutoff parameter (j33|) is field dependent. We therefore follow Ref. [89] and 



find the ground-state binding energy -E^^^^ for the first exciton we are interested in 
here from the transcendental equation 



In 




0. (35) 



In doing so, we first find the exciton Rydberg energy, Ry* (=/ue^/2/i^e^), from this 
equation at F = 0. We use the diameter- and chirality-dependent electron and hole 
effective masses from Ref. pO], and the first bright exciton binding energy of 0.76 eV 
for both (11,0) and (10,0) CN as reported in Ref. [21j from ab initio calculations. We 
obtain Ry* = 4.02 eV and 0.57 eV for the (11,0) tube and (10,0) tube, respectively. 
The difference of about one order of magnitude reflects the fact that these are the 
semiconducting CNs of different types — type-I and type-II, respectively, based on 
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Figure 5: (a) Calculated binding energies of the first bright exciton in the (11,0) and 
(10,0) CNs as functions of the perpendicular electrostatic field applied. Solid lines are 
the numerical solutions to Eq. (j35p . dashed lines are the quadratic approximations 
as given by Eq. (j36p . (b) Field dependence of the effective cutoff Coulomb potential 
(j32p in the (11,0) CN. The dimensionless energy is defined as [Energy]/ 2 jq, according 
to Eq. 

(2n + m) families [90]. The parameters Ry* thus obtained are then used to find 
li^^^^^l as functions of F by numerically solving Eq. ()35p with zo{l,F) given by 
Eqs. (f33D and (l3iD. 

The calculated (negative) binding energies are shown in Fig. O^a) by the solid 
lines. Also shown there by dashed lines are the functions 

Ei''\F) ^ Ei''\l - A,{F)] (36) 

with Ai(F) given by Eq. (l34l) . They are seen to be fairly good analytical (quadratic 
in field) approximations to the numerical solutions of Eq. (j35p in the range of not too 
large fields. The exciton binding energy decreases very rapidly in its absolute value 
as the field increases. Fields of only ~ 0.1 — 0.2 V//um are required to decrease l-E^^^^I 
by a factor of ~ 2 for the CNs considered here. The reason is the perpendicular 
field shifts up the "bottom" of the effective potential (j32]) as shown in Fig. EJ^b) 
for the (11,0) CN. This makes the potential shallower and pushes bound excitonic 
levels up, thereby decreasing the exciton binding energy in its absolute value. As 
this takes place, the shape of the potential does not change, and the longitudinal 
relative electron-hole motion remains finite at all times. As a consequence, no tunnel 
exciton ionization occurs in the perpendicular field, as opposed to the longitudinal 
electrostatic field (Franz-Keldysh) effect studied in Ref. [40j where the non-zero field 
creates the potential barrier separating out the regions of finite and infinite relative 
motion and the exciton becomes ionized as the electron tunnels to infinity. 

The binding energy is only the part of the exciton excitation energy ([7]). Another 
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part comes from the band gap energy ([8]) , where £e and are given by the solutions 
of Eqs. ([27|) and ([28l) . respectively. Solving them to the leading (second) order 
perturbation theory approximation in the field (Appendix D), one obtains 



1 



meAj{F) _ rrihAjiF) 



(37) 



where the electron and hole subband shifts are written separately. This, in view of 
Eq. (|34p . yields the first band gap field dependence in the form 



£;ai)(F)«i^f) 



1 



Ai(F) 



(38) 



The bang gap decrease with the field in Eq. (|38p is stronger than the opposite effect 
in the negative exciton binding energy given (to the same order approximation in 
field) by Eq. pUj) . Thus, the first exciton excitation energy d?]) will be gradually 
decreasing as the perpendicular field increases, shifting the exciton absorption peak 
to the red. This is the basic feature of the quantum confined Stark effect observed 
previously in semiconductor nanomaterials [HH [85l [861 [57] . The field dependences 
of the higher interband transitions exciton excitation energies are suppressed by the 
rapidly (quadratically) increasing azimuthal quantization numbers in the denomi- 
nators of Eqs. (1341) and (j37l). 

Lastly, the perpendicular field dependence of the interband plasmon resonances 
can be obtained from the frequency dependence of the axial surface conductivity 
due to excitons (see Ref. [67J and refs. therein). One has 



-ihwfj 



1,2,... [^1 



exc I ■ 



(39) 



where fj and r are the exciton oscillator strength and relaxation time, respectively. 
The plasmon frequencies are those at which the function Re[l/cT|^(w)] has max- 

(11) (22) 

ima. Testing it for maximum in the domain Eexc <hw < Eexc , one finds the first 
interband plasmon resonance energy to be (in the limit r— )-cxd) 



^(11) ^ / [-E'exc'']^ + [Eexi'] 



.(22) 1 



(40) 



Using the field dependent Eex]} given by Eqs. 0, (l36l) and (l38l) . and neglecting the 

(22) 

field dependence of Eexc , one obtains to the second order approximation in the field 



<)(F)«<) 



(ii)/„p(ii) 



1 +E'g''>/2E, 
1 I ^(22Wii) 

^ ^ J-'exc I ^exc 



Ai(F) 



(41) 



Figure [6] shows the results of our calculations of the field dependences for the 
first bright exciton parameters in the (11,0) and (10,0) CNs. The energy is measured 
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Figure 6: (a),(b) Calculated dependences of the first bright exciton parameters in the 
(11,0) and (10,0) CNs, respectively, on the electrostatic field applied perpendicular to 
the nanotube axis. The dimensionless energy is defined as [Energy] /2^q, according 
to Eq. (|17|) . The energy is measured from the top of the first unperturbed hole 
subband. 



from the top of the first unperturbed hole subband (as shown in Fig. [21 right panel). 
The binding energy field dependence was calculated numerically from Eq. (135p as 
described above [shown in Fig. [5] (a)]. The band gap field dependence and the plas- 
mon energy field dependence were calculated from Eqs. (f37l) and (flTI) . respectively. 
The zero-field excitation energies and zero-field binding energies were taken to be 
those reported in Ref. [51] and in Ref. [21], respectively, and we used the diameter- 
and chirality-dependent electron and hole effective masses from Ref. |90] . As is 
seen in Fig. [6] (a) and (b), the exciton excitation energy and the interband plas- 
mon energy experience red shift in both nanotubes as the field increases. However, 
the excitation energy red shift is very small (barely seen in the figures) due to the 
negative field dependent contribution from the exciton binding energy. So, Eix^c (F) 
and E^^\f) approach each other as the field increases, thereby bringing the total 
exciton energy ([6]) in resonance with the surface plasmon mode due to the non-zero 
longitudinal kinetic energy term at finite temperature]! Thus, the electrostatic field 
applied perpendicular to the CN axis (the quantum confined Stark effect) may be 
used to tune the exciton energy to the nearest interband plasmon resonance, to put 
the exciton-surface plasmon interaction in small-diameter semiconducting CNs to 
the strong-coupling regime. 

^ We are based on the zero-exciton-temperature approximation in here [91] , which is weU justified 
because of the exciton excitation energies much larger than fcsT in CNs. The exciton Hamiltonian 
([4]) does not require the thermal averaging over the exciton degrees of freedom then, yielding the 
temperature independent total exciton energy ([6)|. One has to keep in mind, however, that the 
exciton excitation energy can be affected by the enviromental effect not under consideration in here 
(see Ref. [92]). 
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3.3 The optical absorption 

Here, we analyze the longitudinal exciton absorption line shape as its energy is tuned 
to the nearest interband surface plasmon resonance. Only longitudinal excitons (ex- 
cited by light polarized along the CN axis) couple to the surface plasmon modes 
as discussed at the very beginning of this section (see Ref. [7^ for the perpendic- 
ular light exciton absorption in CNs). We start with the linear (weak) excitation 
regime where only single-exciton states are excited, and follow the optical absorp- 
tion/emission lineshape theory developed recently for atomically doped CNs |10j . 
(Obviously, the absorption line shape coincides with the emission line shape if the 
monochromatic incident light beam is used in the absorption experiment.) Then, the 
non-linear (strong) excitation regime is considered with the photonduced excitation 
of biexciton states. 

When the /-internal state exciton is excited and the nanotube's surface EM 
field subsystem is in vacuum state, the time-dependent wave function of the whole 



Here, |{l/(k)})ea; is the excited single-quantum Fock state with one exciton and 
|{l(k, w)}) is that with one surface photon. The vacuum states are |{0})ex and |{0}) 
for the exciton subsystem and field subsystem, respectively. The coefficients C/(k, t) 
and C(k, cj,t) stand for the population probability amplitudes of the respective 
states of the whole system. The exciton energy is of the form Sj(k) = £'j(k)— ih/r 
with Ef{h) given by Eq. ([6]) and r being the phenomenological exciton relaxation 
time constant [assumed to be such that h/T<^Ef(k)] to account for other possible 
exciton relaxation processes. From the literature we have Tp/j ~ 30 — 100 fs for the 
exciton- phonon scattering [40], ~ 50 ps for the exciton scattering by defects [25l 
[28] . and Trad ~ lOps — 10ns for the radiative decay of excitons |5T]. Thus, the 
scattering by phonons is the most likely exciton relaxation mechanism. 

Using Eqs.([5]) and ([9]), we transform the total Hamiltonian (fT|)- (fTO]l to the k- 
representation (see Appendix A), and apply it to the wave function in Eq. ()12]) . 
We obtain the following set of the two simultaneous differential equations for the 
coefficients Cf(k,t) and C(k,w,t) from the time dependent Schrodinger equation 




(42) 




(43) 



C(k',uj,t)e 



'kk' 



^j;[gW(k,kV)]*C;(k,t)e- 



iEf{k)t/h 



^See the footnote on page 1121 above. 
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The (5-symbol on the left in Eq. (j44p ensures that the momentum conservation is 
fulfilled in the exciton-photon transitions, so that the annihilating exciton creates the 
surface photon with the same momentum and vice versa. In terms of the probability 
amplitudes above, the exciton emission intensity distribution is given by the final 
state probability at very long times corresponding to the complete decay of all 
initially excited excitons, 

/(a;) = |C(k,a;,t^oo)p = lj;|gW(k,k,u;)|2 

/ 



/ 

Jo 



dt'Cf{k,t')e 



-i[Ef{]i)~hw]t'/h 



(44) 



Here, the second equation is obtained by the formal integration of Eq. (I44p over 
time under the initial condition C (k,a;,0) = 0. The emission intensity distribution 
is thus related to the exciton population probability amplitude C/(k, t) to be found 
from Eq. (fl3]l . 

The set of simultaneous equations (03]) and ([H]) [and Eq. (jM]), respectively] 
contains no approximations except the (commonly used) neglect of many-particle 
excitations in the wave function (j42p . We now apply these equations to the exciton- 
surface-plasmon system in small-diameter semiconducting CNs. The interaction 
matrix element in Eqs. (|43p and (|44p is then given by the k-transform of Eq. (jl3p . 
and has the following property (Appendix C) 



^|g(+)(k,k,^)P = l.rl{x)p{x) (45) 



with Tq{x) and p{x) given by Eqs. (fTSj) and (HHD, respectively. We further substitute 
the result of the formal integration of Eq. (jlH) [with C (k, w, 0) = 0] into Eq. ([iS]) , use 
Eq. (USD with p{x) approximated by the Lorentzian calculate the integral over 
frequency analytically, and differentiate the result over time to obtain the following 
second order ordinary differential equation for the exciton probability amplitude 
[dimensionless variables, Eq. (fT7|) ] 

CfifS) + [Axp -Aef + i{xp - Sf)] C0) + {Xf/2fCf{f3) = 0, 

where X/ = [2Axpf /(xp)]^/^ with ff{xp) = rl{xp)p{xp), Asf = h/2-foT, /? = 2jot/fi 
is the dimensionless time, and the k-dependence is omitted for brevity. When the 
total exciton energy is close to a plasmon resonance, £/ ~ Xp, the solution of this 
equation is easily found to be 
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where 6x = Axp — Ae / > and X 



[2AxpTf{ef)]^/'^. This solut ion is vahd when 



Ef^Xp regardless of the strength of the exciton-surface-plasmon couphng. It yields 



If, 



the exponential decay of the excitons into plasmons, \Cf{(3)\ ^ exp[—Tf{ef) 
in the weak coupling regime where the coupling parameter (Xf/Sx)'^ <C 1. 
on the other hand, {Xf/6x)'^ ^ 1, then the strong coupling regime occurs, and 
the decay of the excitons into plasmons proceeds via damped Rabi oscillations, 
|C/(/3)p ~ exp(— cos^(Xj/3/2). This is very similar to what was earlier re- 
ported for an excited two-level atom near the nanotube surface [58l EHl [Ml IS] . Note, 
however, that here we have the exciton-phonon scattering as well, which facilitates 
the strong exciton-plasmon coupling by decreasing 6x in the coupling parameter. In 
other words, the phonon scattering broadens the (longitudinal) exciton momentum 
distribution [93], thus effectively increasing the fraction of the excitons with Sf^Xp. 

In view of Eqs. (I45p and (j46p . the exciton emission intensity (I44p in the vicinity 
of the plasmon resonance takes the following (dimensionless) form 



I(x) 



/•oo 

loisf) V / d(3 Cf{(3) e*(^-^/+^^^/)/5 
Jo 



(47) 



where I{x) = 2^oI{u)/h and Io = T f{ef)/2TT. After some algebra, this results in 



I{x) 



Io{ef)[{x-eff + Axl] 



[{x - EfY - Xj/4]^ + {x- EffiAxl + Ae^ 



2^ ' 



(48) 



where Ax"^ > Ae'j. The summation sign over the exciton internal states is omitted 
since only one internal state contributes to the emission intensity in the vicinity of 
the sharp plasmon resonance. 

The line shape in Eq. (I48p is mainly determined by the coupling parameter 



{X f / Axp)'^ . It is clearly seen to be of a symmetric two-peak structure in the strong 
coupling regime where (Xf/Axp)'^ ^ 1. Testing it for extremum, we obtain the 
peak frequencies to be 



Xi^2 = £/ ± 



Xf 



A 



1 + 



AXr 



X 



f 



AXr 



X 



[terms (Axp)^(Aey)^/Xj are neglected], with the Rabi splitting xi —X2^Xf. In 
the weak coupling regime where {X j / Axp)"^ <^ 1, the frequencies xi and X2 become 
complex, indicating that there are no longer peaks at these frequencies. As this 
takes place, Eq. (jl8|) is approximated with the weak coupling condition, the fact 

2AxpT f{ef), to yield the Lorentzian 

/o(g/)/[l + {Asf/Axpr] 



that xr^Ef, and Xj - 

i{x) 



eff+ ff{ef)/2Jl + {AeJAx, 
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peaked at x = £f, whose half-width-at-half-maximum is slightly narrower, however, 
than Tf{ef)/2 it should be if the exciton-plasmon relaxation were the only relax- 
ation mechanism in the system. The reason is the competing phonon scattering takes 
excitons out of resonance with plasmons, thus decreasing the exciton-plasmon re- 
laxation rate. We therefore conclude that the phonon scattering does not affect the 
exciton emission/absorption line shape when the exciton-plasmon coupling is strong 
(it facilitates the strong coupling regime to occur, however, as was noticed above), 
and it narrows the (Lorentzian) emission/absorption line when the exciton-plasmon 
coupling is weak. 

The non-linear optical susceptibility is proportional to the linear optical response 
function under resonant pumping conditions [94]. This allows us to use Eq. (I48[) to 
investigate the non-linear excitation regime with the photoinduced biexciton for- 
mation as the exciton energy is tuned to the nearest interband plasmon resonance. 
Under these conditions, the third-order longitudinal CN susceptibility takes the 
form ESI M 



I{x) 



1 



1 



e/ + i(r//2 + Aej) 



(ey-kf^|)+.(r//2 + Ae;) 



(49) 

where is the (negative) dimensionless binding energy of the biexciton composed 
of two /-internal state excitons, and xo is the frequency-independent constant. The 
first and second terms in the brackets represent bleaching due to the depopula- 
tion of the ground state and photoinduced absorption due to exciton-to-biexciton 
transitions, respectively. 

The binding energy of the biexciton in a small-diameter (~ 1 nm) CN can be 
evaluated by the method pioneered by Landau [95], Gor'kov and Pitaevski [96] . 
Holstein and Herring |97j — from the analysis of the asymptotic exchange coupling 
by perturbation on the configuration space wave function of the two ground-state 
one-dimensional (ID) excitons. Separating out circumferential and longitudinal de- 
grees of freedom of each of the excitons by means of Eq. (I26p . one arrives at the 
biexciton Hamiltonian of the form [see Fig. [7] (a)] 



1/92 92 \ 



(50) 



+ 



1 



Z2 + AZ| + Zo 

1 



+ 



1 



+ 



1 



\{zi+ Z2)/2 + /\Z\ + Zo 

1 



+ 



Z2\ + Zo \zi - AZ| + Zo 

1 

\{zi+Z2)/2- AZ\+Zo 

1 



\{zi - Z2)/2 + AZ| + zo \{zi - Z2)/2 - AZ\ + zo ' 

Here, z\^2 = -2ei,2 — -2^1,2 is the electron-hole relative motion coordinates of the two 
ID excitons, zo is the cut-off parameter of the effective longitudinal electron-hole 
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Figure 7: (a) Schematic (arbitrary units) of the exchange couphng of two ground- 
state ID excitons to form a biexcitonic state, (b) The couphng occurs in the con- 
figuration space of the two independent longitudinal relative electron-hole motion 
coordinates, zi and Z2, of each of the excitons, due to the tunneling of the system 
through the potential barriers formed by the two single-exciton cusp-type poten- 
tials [bottom, also in (a)], between equivalent states represented by the isolated 
two-exciton wave functions shown on the top. 

Coulomb potential ()32p . and AZ = Z2 — Z\ is the center-of-mass-to-center-of-mass 
inter-exciton separation distance. Equal electron and hole effective masses mg^/i are 
assumed |90j and " atomic units" are used [951 EHl EZ! > whereby distance and energy 
are measured in units of the exciton Bohr radius and in units of the doubled 
ground-state-exciton binding energy lE^, = —2Ry*/vQ, respectively. The first two 
lines in Eq. (j50p represent two isolated non-interacting ID excitons [see Fig. [7] (a)]. 
The last two lines are their exchange Coulomb interactions — electron-hole and 
electron-electron + hole-hole, respectively. 

The Hamiltonian (I50p is effectively two dimensional in the configuration space of 
the two independent relative motion coordinates, zi and 22. Figure [7] (b), bottom, 
shows schematically the potential energy surface of the two closely spaced non- 
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interacting ID excitons [line two in Eq. (jSOp ] in the (-21,^2) space. The surface has 
four symmetrical minima [representing equivalent isolated two-exciton states shown 
in Fig. [7| (b), top], separated by the potential barriers responsible for the tunnel 
exchange coupling between the two-exciton states in the configuration space. The 
coordinate transformation x = {zi — Z2 — AZ)/-v/2, y = {zi + Z2)l\pl places the 
origin of the new coordinate system into the intersection of the two tunnel channels 
between the respective potential minima [Fig. [7](b)], whereby the exchange splitting 
formula of Refs. [951 ESI EZ] takes the form 

C/3,„(AZ)-2Sb = ^J(AZ), (51) 

where Ug^u are the ground and excited state energies, respectively, of the two cou- 
pled excitons (the biexciton) as functions of their center-of-mass-to-center-of-mass 
separation, and 



AZ/V2 

dy 



(52) 

x=0 



is the tunnel exchange coupling integral, where ip{x, y) is the solution to the Schro- 
dinger equation with the Hamiltonian (j50p transformed to the (a;, y) coordinates. 
The factor 2/3! comes from the fact that there are two equivalent tunnel channels 
in the problem, mixing three equivalent indistinguishable two-exciton states in the 
configuration space [one state is given by the two minima on the y-axis, and two 
more are represented by each of the minima on the x-axis — compare Figs. [7] (a) 
and (b)]. 

The function ip{x^y) in Eq. (j52p is sought in the form 

il){x,y) = iljQ{x,y)eyL^[-S{x,y)], (53) 

where -00 = J^o'"'^ 6xp[— (|zi(x, y, AZ)| + |z2(x, y, AZ)|)/fo] is the product of two 
single-exciton wave functionqfj representing the isolated two-exciton state centered 
at the minimum zi = Z2 = (or x = —AZ/\/2, y = 0) of the configuration space 
potential [Fig. [7] (b)], and S{x,y) is a slowly varying function to take into account 
the deviation of tp from i/jo due to the tunnel exchange coupling to another equiv- 
alent isolated two-exciton state centered at zi = AZ, Z2 = —AZ (or x = 
y = 0). Substituting Eq. (j53p into the Schrodinger equation with the Hamiltonian 
(j50p pre-transformed to the (x, y) coordinates, one obtains in the region of interest 

dS f 1 1 1 1 

^0 ^. ^ TTTT^ + 



dx \x + 3AZ/V2 X-AZ/V2 y-V2AZ y + V2AZ 



^This is an approximate solution to the Shrodinger equation with the Hamiltonin given by the 
first two hues in Eq. (|50|) . where the cut-off parameter zq is neglected [89]. This approximation 
greatly simplifies problem solving here, while still remaining adequate as only the long-distance tail 
of tpo is important for the tunnel exchange coupling. 
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up to (negligible) terms of the order of the inter-exciton van der Waals energy and 
up to second derivatives of S. This equation is to be solved with the boundary condi- 
tion S{-AZ/V2,y) = originating from the natural requirement 'ilj{—AZ/^/2,y) = 
i/jo{—AZ/^/2,y), to result in 



S{x,y) = z^o In 



;+3AZ/^/2 



AZ/V2 



2V2AZix+AZ/V2) \ 
+ y^ - 2AZ^ • 



(55) 



After plugging Eqs. (|54p and (|53p into Eq. (|52p . and retaining only the leading 
term of the integral series expansion in powers of uq subject to AZ > 1, Eq. (|5ip 
becomes 

\±^r-Y\ze-'^'/^'> 

The ground state energy Ug of two coupled ID excitons is now seen to go through the 
negative minimum (biexcitonic state) as the inter-exciton center-of-mass-to-center- 
of-mass separation AZ increases (Fig. [8]). The minimum occurs at AZq = 1^0/2, 
whereby the biexciton binding energy is Exx~ (2i^'fe/9z^o)(e/3)^'^°~^, or, expressing 
fo in terms of -Ef, and measuring the energy in units of Ry* , 



i^.vx[.nfl,-l=.-||E.r(|)"^'™ . (56) 

The energy Exx can be affected by the quantum confined Stark effect since \E}\ 
decreases quadratically with the perpendicular electrostatic field applied as shown 
in Fig. O (a). Since e/3 ~ 1, the field dependence in Eq. (|56p mainly comes from 
the pre-exponential factor. So, l-Bxxl will be decreasing quadratically with the field 
as well, for not too strong perpendicular fields. At the same time, the equilibrium 
inter-exciton separation in the biexciton, AZq = vqI2 ~ \E\\~'^I'^ , will be slowly 
increasing with the field consistently with the lowering of |i?jo:|- In the zero field, 
one has roughly Exx ~ l-E'bP''^ ~-^CJV^ ^"-"^ '^'^ biexciton binding energy versus the 
CN radius Rcn ^^h\'^ R'cn' ^ reported in Ref. [19] from variational calculations), 
pretty consistent with the dependence obtained numerically [32j. Interestingly, 
as Rcn goes down, l-Exxl goes up faster than \E\,\ does. This is partly due to the 
fact that AZq slowly decreases as Rcn goes down, — a theoretical argument in 
support of experimental evidence for increased exciton-exciton annihilation in small 



diameter CNs [981 l99l finn] . 

Figure [8] shows the ground state energy Ug{AZ) of the coupled pair of the first 
bright excitons, calculated from Eq. (I55p for the semiconducting (11,0) CN exposed 
to different perpendicular electrostatic fields. The inset shows the field dependences 
of Exx [as given by Eq. (i56]l ] and of AZq. All the curves are calculated using the 
field dependence of Ej, obtained as described in the previous subsection (Figs. [5] 
and [6]). They exhibit typical behaviors discussed above. 

Figure [9] compares the linear response lineshape (|48p with the imaginary part 
of Eq. ()49[) representing the non-linear optical response function under resonant 
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Figure 8: Calculated ground state energy Ug of the coupled pair of the first bright 
excitons in the (11,0) CN as a function of the center-of-mass-to-center-of-mass inter- 
exciton distance AZ and perpendicular electrostatic field applied. Inset shows 
the biexciton binding energy Exx and inter-exciton separation AZq (y- and x- 
coordinates, respectively, of the minima in the main figure) as functions of the field. 



pumping, both calculated for the 1st bright exciton in the (11,0) CN as its energy 
is tuned (by means of the quantum confined Stark effect) to the nearest plasmon 
resonance (vertical dashed line in the figure). The biexciton binding energy in 
Eq. ([I9|) was taken to be Exx ~ 52 meV as given by Eq. (j56|) in the zero field. 
[Weak field dependence of Exx (inset in Fig. [8]) plays no essential role here as 
l-E'xxl \Eb\ ~ 0.76 eV regardless of the field strength.] The phonon relaxation 
time Tph = 30 fs was used as reported in Ref. [29], since this is the shortest one 
out of possible exciton relaxation processes, including exciton-exciton annihilation 
(ree~l ps tSSj). Clear line (Rabi) splitting effect ~0.1 eV is seen both in the linear 
and in non-linear excitation regime, indicating the strong exciton-plasmon coupling 
both in the single-exciton states and in the biexciton states as the exciton energy is 
tuned to the interband surface plasmon resonance. The splitting is not masked by 
the exciton-phonon scattering. 

This efi^ect can be used for the development of new tunable optoelectronic de- 
vice applications of optically excited small-diameter semiconducting CNs in areas 
such as nanophotonics, nanoplasmonics, and cavity quantum electrodynamics, in- 
cluding the strong excitation regime with optical non-linearities. In the latter case. 
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Figure 9: [(a), (b), and (c)] Linear (top) and non- linear (bottom) response functions 
as given by Eq. (|18|) and by the imaginary part of Eq. (09]), respectively, for the 
first bright exciton in the (11,0) CN as the exciton energy is tuned to the nearest 
interband plasmon resonance (vertical dashed line). Vertical lines marked as X 
and XX show the exciton energy and biexciton binding energy, respectively. The 
dimensionless energy is defined as [Energy] /2^q, according to Eq. (fTTI) . 

the experimental observation of the non-linear absorption line splitting predicted 
here would help identify the presence and study the properties of biexcitonic states 
(including biexcitons formed by excitons of different subbands [33]) in individual 
single- walled CNs, due to the fact that when tuned close to a plasmon resonance 
the exciton relaxes into plasmons at a rate much greater than T~j^ ( ^t^J^), totally 
ruling out the role of the competing exciton-exciton annihilation process. 



4 Casimir Interaction in Double- Wall 
Carbon Nanotubes 

Here, we consider the Casimir interaction between two concentric cylindrical gra- 
phene sheets comprising a double-wall CN, using the macroscopic QED approach 
employed above to study the exciton-surface-plasmon interactions in single wall 
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nanotubesO The method is fuhy adequate in this case as the Casimir force is 
known to originate from quantum EM field fluctuations. The fundamental nature 
of this force has been studied for many years since the prediction of the attraction 
force between two neutral metallic plates in vacuum (see, Refs. \^\ I56j). After the 
first report of observation of this spectacular effect [lOlj . new measurements with 
improved accuracy have been done involving different geometries [1021 11031 I104j - 
The Casimir force has also been considered theoretically with methods primarily 
based on the zero-point summation approach and Lifshitz theory |105t 1106] . 

The Casimir effect has acquired a much broader impact recently due to its im- 
portance for nanostructured materials, including graphite and graphitic nanostruc- 
tures [56J which can exist in different geometries and with various unique electronic 
properties. Moreover, the efficient development and operation of modern micro- and 
nano-electromechanical devices are limited due to effects such as stiction, friction, 
and adhesion, originating from or closely related to the Casimir effect |57j . 

The mechanisms governing the CN interactions still remain elusive. It is known 
that the system geometry [107^ 1108] and dielectric response f45^, '62] have a profound 
effect on the interaction, in general, but their specific functionalities have not been 
qualitatively and quantitatively understood. Since CNs of virtually the same radial 
size can possess different electronic properties, investigating their Casimir interac- 
tions presents a unique opportunity to obtain insight into specific dielectric response 
features affecting the Casimir force between metallic and semiconducting cylindri- 
cal surfaces. This can also unveil the role of collective surface excitations in the 
energetic stability of multi-wall CNs of various chiral combinations. 

Since Lifshitz theory cannot be easily applied to geometries other than parallel 
plates, researchers have used the Proximity Force Approximation (PFA) to calculate 
the Casimir interaction between CNs [1071 I109j (see also Ref. |56j for the latest 
review). The method is based on approximating the curved surfaces at very close 
distances by a series of parallel plates and summing their energies using the Lifshitz 
result. Thus, the PFA is inherently an additive approach, applicable to objects at 
very close separations (still to be greater than objects inter-atomic distances) under 
the assumption that the CN dielectric response is the same as the one for the plates. 
This last assumption is very questionable as the quasi-lD character of the electronic 
motion in CNTs is known to be of principal importance for the correct description 
of their electronic and optical properties [H [59l IM] . 

We model the double-wall CN by two infinitely long, infinitely thin, continu- 
ous concentric cylinders with radii -Ri,2i immersed in vacuum. Each cylinder is 
characterized by the complex dynamic axial dielectric function ezziRi,2,^) with the 
z-direction along the CN axis as shown in Fig. [101 The azimuthal and radial com- 
ponents of the complete CN dielectric tensor are neglected as they are known to be 
much less than Czz for rnost CNs [6l] • The QED quantization scheme in the presence 



In this Section only, the International System of units is used to make the comparison easier of 
our theory with other authors' results. 
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Figure 10: Schematic of the two concentric CNs in vacuum. The CN radii are Ri 
and i?2- The regions between the CN surfaces are denoted as (1), (2), and (3). 



of CNs [491 [62] generates the second-quantized Hamiltonian 

of the vacuum-type medium assisted EM field, with the bosonic operators P and 
/ creating and annihilating, respectively, surface EM excitations of frequency oj at 
points Ri,2 = {-Ri, 2, ¥'1,2, 2^1,2} of the double-wall CN system. The Fourier-domain 
electric field operator at an arbitrary point r = (r, (/?, z) is given by 

£'(r, w) = iw/io / (iRiG(r, Rj, w) • J(Ri, w), 

j=l,2 

where G(r, Rj,a;) is the dyadic EM field Green's function (GF), and 

is the surface current density operator selected in such a way as to ensure the 
correct QED equal-time commutation relations for the electric and magnetic field 
operators [l9l[62]. Here, is the unit vector along the CN axis, Eq, hq, and c are the 
dielectric constant, magnetic permeability, and vacuum speed of light, respectively. 
The dyadic GF satisfies the wave equation 

^2 

V X V X G(r, r', lo) - ^ G(r, r', 00) = 6(r - r') I (57) 
with I being the unit tensor. The GF can further be decomposed as follows 



G(-/) = G(0)5., + Gi:£ 
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where G^*^-* and G^^'/^^^ represent the contributions of the direct and scattered waves, 
respectively llllj . with a point-hke field source located in region s and the 
field registered in region / (see Fig. [TOj) . The boundary conditions for Eq. (|57p 
are obtained from those for the electric and magnetic field components on the CN 
surfaces [l5l |59] , which result in 



Gt- X 



G(r,r',a;)|^+^- G(r,r',a;)| 



0, 



(58) 



Br- X V X 



G{r,r',uj)\ +- G(r,r',a;)| 

'-"■1,2 '-"-1,2 



(59) 



where is the unit vector along the radial direction. The discontinuity in Eq. (j59p 
results from the full account of the finite absorption and dispersion for both CNs by 
means of their conductivity tensors cr^^'^^ approximated by their largest components 



-|^[4l'^)(i?,„a;)-l] 



(60) 



[compare with Eq. ([21]) ]. 



Following the procedure described in Refs. [11011111] . we expand G^'^) and G^*'''^'* 



scatt 



into series of even and odd vector cylindrical functions with unknown coefficients to 
be found from Eqs. ([58p and ([59p . This splits the EM modes in the system into TE 
and TM polarizations, with Eqs. ([58p and ([59p yielding a set of 32 equations (16 
for each polarization) with 32 unknown coefficients. The unknown coefficients are 
found determining the dyadic GF in each regionlll 

Using the expressions for the electric and magnetic fields, the electromagnetic 
stress tensor is constructed [l9l 1112] 



T(r, r') = Ti(r, r') + T2(r, r') - -ITr [Ti(r, r') + T2(r, r')] 



(61) 



fc rco 2 

Ti(r,r') = - / da;— Im[G(r,r', 



T2(r,r') 



- — / duj Im 
71" Jo 



V X G(r,r ,L^)x V' 



(62) 



(63) 



We are interested in the radial component T^-r which describes the radiation pressure 
of the virtual EM field on each CN surface in the system. The Casimir force per 
unit area exerted on the surfaces is then given by 



Fi = lim < lim 

r-^Ri I r'^r 



Tl^{ry)-T^'^(r,v') 



i = 1,2 



(64) 



^Due to the lengthy and tedious algebra, this derivation will be presented in a separate longer 
communication. 
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The forces Fi^2 calculated from Eq. are of equal magnitude and opposite di- 
rection, indicating the attraction between the cylindrical surfaces. The Casimir force 
thus obtained accounts simultaneously for the geometrical curvature effects (through 
the GF tensor) and the finite absorption and dissipation of each CN [through their 
dielectric response functions (f60]l ]. The dielectric response functions of particu- 
lar CNs were calculated from the CN realistic band structure as described above, 
in Section [3l We decomposed them into the Drude contribution and the contri- 
bution originating from (transversely quantized) interband electronic transitions, 
^zz = + ^zz*'^'^) ill order to be able to see how much each individual contribution 
affects the inter-tube Casimir attraction. 

It is interesting to consider the case of infinitely conducting parallel plates first 
using Eq. ([Ml) . This is obtaned by taking the limits alz — ^ oo and Ri^2 — ^ oo 
while keeping constant the inter-tube distance, Ri — R2 = d. We find 



X { [xlK^xi) + (n^ + xl) K^ixi)] [ll{xi)Kn{x2)/Kn{xi) - 2/„(xi)/„(x2)] 
- [xll'^ixi) + (n^ + xl) Ilixi)] Kn{xi)Kn{x2) 

- 2 [x?/;(xiX(xi) + (n2 + x\) In{xi)Kn{xi)\ In{x2)Kn{xi)] 



where xi^2 = xRi^2-, In{x) and Kn{x) are the modified Bessel functions of the first 
and second kind, respectively. The above expression is obtained by making the tran- 
sition to imaginary frequencies a; — t- io;, and using the Euclidean rotation technique 
as described in Refs. |1121 1113j . This can further be evaluated by summing up the 
series over n using the large-order Bessel function expansions |114j . This results 
in F ~ (— l/3)(/ic7r^/240(i'') which is 1/3 of the well-known result for two paral- 
lel plates |491 [56] . This deviation originates from ezz 7^ only and the remaining 
dielectric tensor components being zero in our model. 

Figure [TT] presents results from the numerical calculations of F as a function of 
the inter-tube surface-to-surface distance for various pairs of CNs with their realistic 
chirality dependent dielectric responses taken into account. We have chosen the 
inner CN to be the achiral (12, 12) metallic nanotube, and to change the outer 
tubes. As R2 is varied, one can envision double wall CNs consisting of metal/metal 
or metal/semiconductor combinations of different chiralities but of similar radial 
dimensions. 

Figure [11] shows that F decreases in strength as the surface-to-surface distance 
increases. This dependence is monotonic for the zigzag (m, 0) and armchair (n, n) 
outer tubes, but it happens at different rates. The attraction is stronger if the outer 
CN is an armchair (n, n) one as compared to the attraction for the outer (m, 0) 
nanotubes. At the same time, for chiral tubes the Casimir force decreases as a 
function of d in a rather irregular fashion. It is seen that for relatively small d. 
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Figure 11: The Casimir force per unit area as a function of the inter-tube separation 
d, for different pairs of CNs. The inset shows force found with the full dielectric 
function and the Drude contribution only for the same CN pairs indicated in the 
figure. 



the interaction force can be quite different. For example, the attraction between 
(27,4)@(12,12) and (21, 13)@(12, 12) differ by ~ 20 % in favor of the second pair, 
even though the radial difference is only 0.2 A. The differences between the different 
CNs become smaller as their separation becomes larger, and they eventually become 
negligible as the Casimir force diminishes at large distances. 

We also calculate the Casimir force using the e^(w) contribution alone in each 
dielectric function. The inset in Fig. [TT] indicates that the attraction is stronger 
when the interband transitions are neglected. The decay of -F as a function of 
d is monotonic. Including the e^"**^^(a;) term not only reduces the force, but also 
introduces non-linearities due to the chirality dependent optical excitations. At large 
surface-to-surface separations, the discrepancies between the force calculated with 
the full dielectric response, and those obtained with the Drude term only become 
less significant. We find that for (i~15 A, this difference is less than 10 %. 

To investigate further the important functionalities originating from the cylindri- 
cal geometry and the CN dielectric response properties, F is calculated for different 
achiral inner/outer nanotube pairs. Studying zigzag and armchair CNs allows track- 
ing generalities from e{uj) in a more controlled manner. The results are presented in 
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Figure 12: The Casimir force per unit area as a function of the inter-tube separation 
d for selected CN pairs. The insets show the EELS spectra for several CNs. 

Fig. 1121 We have chosen representatives of three inner CN types - metallic (12, 12), 
semi-metallic (21,0), and semiconducting (20,0) tubules. They are of similar radii, 
8.14 A, 8.22 A, and 7.83 A, respectively. We see that depending on the outer nan- 
otube types, the F versus d curves are positioned in three groups. The weakest 
interaction is found when there are two zigzag concentric CNs (top two curves). 
The fact that some of these are semi-metallic and others are semiconducting does 
not seem to influence the magnitude and monotonic decrease of the Casimir force. 

The attraction is stronger when there is a combination of an armchair and a 
zigzag CNT as compared to the previous case. The curves for (m, 0)@(12, 12), 
(n, n)@(21, 0), and (n, n)@(20, 0) are practically overlapping, meaning that the spe- 
cific location of the zigzag and armchair tubes (inner or outer) is of no significance 
to the force. The small deviations can be attributed to the small differences in the 
inner CN radii. Finally, we see that the strongest interaction occurs between two 
armchair CNs (red curve). These functionalities are not unique just for the con- 
sidered CNs. We have performed the same calculations for many different achiral 
tubes, and we always find that the strongest interaction occurs between two arm- 
chair CNs and the weakest — between two zigzag CNs (provided that their radial 
dimensions are similar). 
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The results from these calculations are strongly suggestive that the CN collec- 
tive excitation properties have a strong effect on their mutual interaction. This is 
particularly true for the relatively small distances of interest here, for which the 
dominant contribution of plasmonic modes to the Casimir interactions has been 
realized for planar |115j and linear |117j metallic systems. To elucidate this issue 
here, we calculate the EELS spectra, given by — Im[l/e(a;)], and compare them for 
various inner and outer CNs combinations — Fig. [12] (inset). 

Considering F as a function of d and the specific form of the EELS spectra, it 
becomes clear from the inset in Fig. [12] that the low frequency plasmon excitations, 
given by peaks in — Im[l/e(Li;)], are key to the strength of the Casimir force. We 
always find that the strongest force is between the tubules with well pronounced over- 
lapping low frequency plasmon excitations. This is consistent with the conclusion 
of Ref. |117j for generic ID-plasmonic structures. However, in our case we deal with 
the interband plasmons originating from the space quantization of the transverse 
electronic motion, and, therefore, having quite a different frequency-momentum 
dispersion law (constant) as compared to that normally assumed (linear) for plas- 
mons [50j. A weaker force is obtained if only one of the CNs supports strong low 
frequency interband plasmon modes. The weakest interaction happens when neither 
CN has strong low frequency plasmons. For the cases shown in Fig [T2l one finds 
well pronounced overlapping plasmon transitions in the (12, 12) CN at uji = 2.18 eV 
and UJ2 = 3.27 eV, and at ui = 1.63 eV and U2 = 2.45 eV in the (17, 17) CN. At 
the same time, no such well defined strong low frequency excitations in the (21,0) 
and (30,0) CNs are found. Figure [l2] shows that the attraction in (17, 17)@(12, 12) 
is much stronger than the attraction in (30, 0)@(21, 0), even though the radial sizes 
of the involved CNs are approximately the same. One also notes that for the case 
of (17, 17)@(21,0) there is only one such low frequency excitation coming from the 
armchair tube and, consequently, the Casimir force has an intermediate value as 
compared to the above discussed two cases. 

We performed calculations of the Casimir force between many CN pairs and 
made comparisons between the relevant regions of the EELS spectra. It is found 
that, in general, armchair tubes always have strong, well pronounced interband 
plasmon excitations in the low frequency range. Zigzag and most chiral CNs have 
low frequency interband plasmons [37], too, but they are not as near as well pro- 
nounced as those in armchair tubes; their stronger plasmon modes are found at 
higher frequencies. 

These studies are indicative of the significance of the collective response prop- 
erties of the involved CNs. Specifically, the collective low energy plasmon excita- 
tions and their relative location can result in nanotube attraction with different 
strengths. We further investigate this point by considering a double wall CN with 
radii Ri = 11.63 Aand R2 = 8.22 A. The dielectric function of each tube is taken 
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Figure 13: The Casimir force per unit area as a function of the outer CN plasmon 
frequency, while the inner CN plasmon peak uj2 is constant. Results are shown for 
four values of u}2- The dielectric functions are modeled by a generic Lorentzian as 
given by Eq. (f65]l . 



to be of the generic Lorentzian form 

ezz{Ri,2,uj) = 1 2 2 — r~~F 

with the typical for nanotubes values = 2.7 eV and T = 0.03 eV [l5]. Then, the 
EELS spectrum has only one plasmon resonance at u}i^2 for each tube. This generic 
form allows us to change the relative position and strength of the plasmon peaks 
and uncover more characteristic features originating from the EELS spectra. 

In Fig. \T3\ the force as a function of plasmon frequency resonances of the outer 
CN is shown when the plasmon transition for the inner CNT is kept constant (four 
values are chosen for u}2)- One sees that the local minima in F versus u) occur 
when ui and U2 coincide. In fact, the strongest attraction happens when both CNs 
have the lowest plasmon excitations at the same frequency uji = UJ2 = 0.81 eV. It 
is evident that the existence of relatively strong low frequency EELS spectrum and 
an overlap between the relevant plasmon peaks of the two structures is necessary to 
achieve a strong interaction. 
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This study clearly demonstrates the crucial importance of the collective low en- 
ergy surface plasmon excitations at relatively close surface-to-surface separations 
along with the cylindrical circular geometry of the double-wall CN system. The 
QED approach we used provides the unique opportunity to investigate these fea- 
tures together, or separately, and to uncover underlying mechanisms of the energetic 
stability of different double- wall CN combinations. An additional advantage here is 
that we can calculate the dielectric function explicitly for each chirality. Thus, we 
can determine unambiguously how the semiconducting or metallic nature of each 
CN contributes to their mutual interaction. 

5 Conclusion 

We have shown that the strong exciton-surface-plasmon coupling effect with char- 
acteristic exciton absorption line (Rabi) splitting ~0.1 eV exists in small-diameter 
(< 1 nm) semiconducting CNs. The splitting is almost as large as the typical exci- 
ton binding energies in such CNs (~0.3 — 0.8 eV [IHl EOj [Ml EI]), and of the same 
order of magnitude as the exciton-plasmon Rabi splitting in organic semiconductors 
(~ 180 meV [53]). It is much larger than the exciton-polariton Rabi splitting in 
semiconductor microcavities (~ 140 — 400 fieV [HTl |82l |83] ) , or the exciton-plasmon 
Rabi splitting in hybrid semiconductor- metal nanoparticle molecules [54j . 

Since the formation of the strongly coupled mixed exciton-plasmon excitations 
is only possible if the exciton total energy is in resonance with the energy of an 
interband surface plasmon mode, we have analyzed possible ways to tune the ex- 
citon energy to the nearest surface plasmon resonance. Specifically, the exciton 
energy may be tuned to the nearest plasmon resonance in ways used for the ex- 
citons in semiconductor quantum microcavities — thermally (by elevating sample 
temperature) (HH [82l [83], and/or electrostatically [H [85l [Ml [87] (via the quan- 
tum confined Stark effect with an external electrostatic field applied perpendicular 
to the CN axis). The two possibilities influence the different degrees of freedom of 
the quasi- ID exciton — the (longitudinal) kinetic energy and the excitation energy, 
respectively. 

We have studied how the perpendicular electrostatic field affects the exciton 
excitation energy and interband plasmon resonance energy (the quantum confined 
Stark effect). Both of them are shown to shift to the red due to the decrease in 
the CN band gap as the field increases. However, the exciton red shift is much 
less than the plasmon one because of the decrease in the absolute value of the 
negative binding energy, which contributes largely to the exciton excitation energy. 
The exciton excitation energy and interband plasmon energy approach as the field 
increases, thereby bringing the total exciton energy in resonance with the plasmon 
mode due to the non-zero longitudinal kinetic energy term at finite temperature. 

The noteworthy point is that the strong exciton-surface-plasmon coupling we 
predict here occurs in an individual CN as opposed to various artificially fabricated 
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hybrid plasmonic nanostructures mentioned above. We strongly believe this phe- 
nomenon, along with its tunability feature via the quantum confined Stark effect we 
have demonstrated, opens up new paths for the development of CN based tunable 
optoelectronic device applications in areas such as nanophotonics, nanoplasmonics, 
and cavity QED. One straightforward application like this is the CN photolumi- 
nescence control by means of the exciton-plasmon coupling tuned electrostatically 
via the quantum confined Stark effect. This complements the microcavity controlled 
CN infrared emitter application reported recently [27J, offering the advantage of less 
stringent fabrication requirements at the same time since the planar photonic micro- 
cavity is no longer required. Electrostatically controlled coupling of two spatially 
separated (weakly localized) excitons to the same nanotube's plasmon resonance 
would result in their entanglement |1H \T2\ [T3] , the phenomenon that paves the way 
for CN based solid-state quantum information applications. Moreover, CNs com- 
bine advantages such as electrical conductivity, chemical stability, and high surface 
area that make them excellent potential candidates for a variety of more practical 
applications, including efficient solar energy conversion energy storage yL4j, and 
optical nanobiosensorics 02]. However, the photoluminescence quantum yield of 
individual CNs is relatively low, and this hinders their uses in the aforementioned 
applications. CN bundles and films are proposed to be used to surpass the poor per- 
formance of individual tubes. The theory of the exciton-plasmon coupling we have 
developed here, being extended to include the inter-tube interaction, complements 
currently available 'weak-coupling' theories of the exciton-plasmon interactions in 
low-dimensional nanostructures |54i I121j with the very important case of the strong 
coupling regime. Such an extended theory (subject of our future publication) will 
lay the foundation for understanding inter-tube energy transfer mechanisms that 
affect the efficiency of optoelectronic devices made of CN bundles and films, as well 
as it will shed more light on the recent photoluminescence experiments with CN 
bundles |43l [13] and multi- walled CNs [122j , revealing their potentialities for the 
development of high- yield, high-performance optoelectronics applications with CNs. 

In addition, we have first applied the macroscopic QED approach suitable for 
dispersing and absorbing media to study the Casimir interaction in a double-wall 
carbon nanotube systems with the realistic dielectric response taken into account. 
We found that at distances similar to the equilibrium separations between graphitic 
surfaces (~ 3 A) , the attraction is dominated by the low energy (interband) plasmon 
excitations of both CNs. The key attributes of the EELS spectra are the existence of 
low frequency plasmons, their strong and well pronounced nature, and the overlap 
between the low frequency plasmon peaks belonging to the two CNTs. Thus, the 
chiralities of concentric graphene sheets with similar radial sizes exhibiting these fea- 
tures will be responsible for forming the most preferred CN pairs. As the inter-tube 
separation increases, the plasmon effect diminishes and the collective excitations 
originating from the nanotube metallic or semiconducting nature do not influence 
the interaction in a profound way. 
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We expect our results to pave the way for the development of new generation of 
tunable optoelectronic and nano-electromechanical device applications with single- 
wall and multi-wall carbon nanotubes. 
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Appendix A 

Exciton interaction with the surface EM field 

We follow our recently developed QED formalism to describe vacuum-type EM 
effects in the presence of quasi-lD absorbing and dispersive bodies [SU \59\ EOl EH 
l62t [9]. The treatment begins with the most general EM interaction of the surface 
charge fluctuations with the quantized surface EM field of a single-walled CN. No 
external field is assumed to be applied. The CN is modelled by a neutral, infinitely 
long, infinitely thin, anisotropically conducting cylinder. Only the axial conductivity 
of the CN, (T22, is taken into account, whereas the azimuthal one, (7^^,^, is neglected 
being strongly suppressed by the transverse depolarization effect [631 EH ES EH EZl 
I68j . Since the problem has the cylindrical symmetry, the orthonormal cylindrical 
basis {er,e^,ez} is used with the vector e^ directed along the nanotube axis as 
shown in Fig. [H The interaction has the following form (Gaussian system of units) 



-A(n + fL'^: 



Hint = H^t + H^t (66) 



pi:^-gA(n + ri;OJ +Y.q^^{^ + v\ 



rriiC 

n,i n,i 



where c is the speed of light, m^, qi, !•„ , and pn are, respectively, the masses, 
charges, coordinate operators and momenta operators of the particles (electrons 
and nucleus) residing at the lattice site n = Iln = {RcNjV'n, Zn} associated with a 
carbon atom (see Fig. [T]) on the surface of the CN of radius Rc n ■ The summation 
is taken over the lattice sites, and may be substituted with the integration over the 
CN surface using Eq. The vector potential operator A and the scalar potential 
operator represent the nanotube's transversely polarized and longitudinally po- 
larized surface EM modes, respectively. They are written in the Schrodinger picture 
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as follows 

A(n) = / doj—E'^{n,uj) + h.c., (67) 



/■°° " II 

Vn(^(n)= j duj'E (n, w) + /i.e.. 



We use the Coulomb gauge whereby Vn • A(n) = 0, or [pn\ A(n + rn^)] = 0. 

The total electric field operator of the CN-modified EM field is given for an 
arbitrary r in the Schrodinger picture by 

E(r) = / dujtj^r,uj) +h.c. = [ do; [E"^(r, w) + e" (r, w)] + /i.e. 
Jo Jo 

with the transversely (longitudinally) polarized Fourier-domain field components 
defined as 



E^^"\r,w)= y"dr' 5^(11) (r-r')-E(r», (70) 



where 



'5l,(r) = -V.V,^, (71) 

are the longitudinal and transverse dyadic 5-functions, respectively. The total field 
operator (j69]) satisfies the set of the Fourier-domain Maxwell equations 

V x'E{r,uj) =ikH{r,uj), (72) 

Att - 

V xU(r,u}) = -ikE(r,uj) + — I(r,w), (73) 

c 

where H = {ik)^^'V x E is the magnetic field operator, k = uj/c, and 

i(r,a;) = ^5(r-n)J(n,a;), (74) 

n 

is the exterior current operator with the current density defined as follows 



J(n,u;) = W /(n,a;)e2 (75) 

to ensure preservation of the fundamental QED equal-time commutation relations 
(see, e.g., ^Tj) fo'^ the EM field components in the presence of a CN. Here, azz is the 
CN surface axial conductivity per unit length, and /(n,a;) along with its counter- 
part /^(n, w) are the scalar bosonic field operators which annihilate and create. 
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respectively, single-quantum EM field excitations of frequency u at the lattice site 
n of the CN surface. They satisfy the standard bosonic commutation relations 

[/(n,^),/"^(m,w')] = ^nmS{(^ -(^'), (76) 

[/(n,a;),/(m,a;')] = [f (n, c^), /t(m, c^')] = 0. 
One further obtains from Eqs. (|72| ) -(|75 p that 

E(r,a;) = iA;^^G(r,n,a;)-J(n,u;), (77) 



and, according to Eqs. ([69|) and ([70 
-(II), X .,47r 



E^"^(r,a;) = iyt— ^^(ll)G(r,n,cj)-J(n,w), (78) 



c 

n 



where """G and "G are the transverse part and the longitudinal part, respectively, 
of the total Green tensor G = """G + "G of the classical EM field in the presence of 
the CN. This tensor satisfies the equation 

(V X V X - e)^^ Gaz (r, n,oj) = 6{r- n) (79) 

a=r,ip,z 

together with the radiation conditions at infinity and the boundary conditions on 
the CN surface. 

All the 'discrete' quantities in Eqs. (f7il) - (f79]l may be equivalently rewritten in 
continuous variables in view of Eq. ([3]). Being applied to the identity 1 = "^^Snm, 
Eq. da]) yields 

5nm = Sod{Kn-Km). (80) 

This requires to redefine 

/(n,a;) = v^/(R„,c^), P{ii,u;) = v^/t(R„,a;) (81) 



in the commutation relations (j76p . Similarly, from Eq. (j77p . in view of Eqs. ([3]), (j75p 
and ([5T]) . one obtains 

G(r,n,a;) = v^G(r,R„,,w), (82) 
which is also valid for the transverse and longitudinal Green tensors in Eq. 



Next, we make the series expansions of the interactions H^^^ and H^^^ in Eq. ([66]) 
about the lattice site n to the first non-vanishing terms, 

Mil ^-Y — A(n) • p« + y -^A\n), (83) 

n,i n,i 

« E '^^^n • ^^n\ (84) 



40 



I.V.Bondarev, L.M.Woods, and A.Popescu 



and introduce the single-lattice-site Hamiltonian 



H^ = eo\0){0\+^{eo + hjf)\f){f\ (85) 
/ 

with the completeness relation 

|0)(0|+^|/)(/|=i. (86) 
/ 

Here, Sq is the energy of the ground state |0) (no exciton excited) of the carbon atom 
associated with the lattice site n, eo + hu}f is the energy of the excited carbon atom 
in the quantum state |/) with one /-internal-state exciton formed of the energy 
Eexc = h^f- In view of Eqs. (j85p and ()86p . one has 

.(i) drW rrii (j) rrii ~ (i) £^ f 

at in in 

^^Y.^f{i^\^^^\f)Br.j- {f\itMBij) (87) 
/ 



and 



E((0|rnV)i?n,/ + (/|r[;V)</), (88) 
/ 



where (0|rn^|/) = (/|rn''|0) in view of the hermitian and real character of the coor- 
dinate operator. The operators i3n,/ = 1 0)(/ 1 and i?^ j = |/)(0| create and annihilate, 
respectively, the /-internal-state exciton at the lattice site n, and exciton-to-exciton 
transitions are neglected. In addition, we also have 

5ii<5a/3 = ;^[(p^:^)a,(i-^.''^)/3], (89) 



where a, /3 = r, ip, z. Substituting these into Eqs. ([55]) and dSU [commutator 
goes into the second term of Eq. (j83p which is to be pre-transformed as follows 
j ^ pqiqjA{n)aA{rY)p5ij5aj3/'2'rniC^], one arrives at the following (electric dipole) 
approximation of Eq. (j66p 



Hint - h'^^i + ^i^l (90) 
E^d/.A(n) 



^I.f-^n,/ + 7^d^A(n) 

n,/ 



he 



n,/ 
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with d;^ = (0|dn|/) = (/|dn|0), where dn = YliQi^n^ total electric dipole 

moment operator of the particles residing at the lattice site n. 

The Hamiltonian (j90p is seen to describe the vacuum-type exciton interaction 
with the surface EM field (created by the charge fluctuations on the nanotube sur- 
face). The last term in the square brackets does not depend on the exciton operators, 
and therefore results in the constant energy shift which can be safely neglected. We 
then arrive, after using Eqs. (j67|) . ([68]) . ([75]) . and ([78|) . at the following second 
quantized interaction Hamiltonian 

duj [g'j^\n,m,uj)Blj - g-f\n,m,u})Bnj]f{ui,uj) + h.c, (91) 

n,m,/ - 

where ^ 

gf\n,m,uj) = g|(n,m,w) ± — g''(n,m,a;) (92) 

with 

gjr(")(n,m,L^) = -i^^^hwKea,,{RcN.u:) ^ (d^), ^(ll)C(n, m, w), (93) 



a=r,ifi,z 



and 



-(\\)a 



--ldr6tlj^\n-r)G^,{r,in,co). (94) 



This yields Eqs. (ll0p ~( fT2]) after the strong transverse depolarization effect in CNs 
is taken into account whereby d^ ~ {d(^)zez. 

Appendix B 

Green tensor of the surface EM field 

Within the model of an infinitely thin, infinitely long, anisotropically conducting 
cylinder we utilize here, the classical EM field Green tensor is found by expanding 
the solution to the Green equation ()79p in series in cylindrical coordinates, and 
then imposing the appropriately chosen boundary conditions on the CN surface to 
determine the Wronskian normalization constant (see, e.g., Ref. [118j). 

After the EM field is divided into the transversely and longitudinally polarized 
components according to Eqs. ([69|) - (j7T]) . the Green equation ([79]) takes the form 

^ (Vx Vx -/c^)^^ ^^G„^(r,n,w) +llGa^(r,n,w)] = 5{r - n) (95) 

with the two additional constraints, 

X] V"^G„,(r,n,w) = (96) 
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and 

^ ea/37V/3l'G^2(r,n,a;) = 0, (97) 



where eQ,^^ is the totahy antisymmetric unit tensor of rank 3. Equations ()96p and 
(j97p originate from the divergence-less character (Coulomb gauge) of the transverse 
EM component and the curl- less character of the longitudinal EM component, re- 
spectively. The transverse '^Gaz and longitudinal "Gq,2 Green tensor components 
are defined by Eq. (j94p which is the corollary of Eq. (j70p using the Eqs. (I77p and 
(jTHj) . Equation ([95]) is further rewritten in view of Eqs. ([96j) and ([97|) . to give the 
following two independent equations for '^Gzz and "Gz^ we need 

{A + k^)^GUr,n,oj) = -6iir-n), (98) 
A:2|lG,,(r,n,a;) = -4(r-n) (99) 

with the transverse and longitudinal delta- functions defined by Eq. ()7ip . 

We use the differential representations for the transverse '^Gzz and longitudinal 
^^Gzz Green functions of the following form [consistent with Eq. (j94p ] 



^Gzz{r,n,oo) = (^^V,V^ + 1^ ^(r, n, a;), (100) 

^^Gzzir,n,uj) = -^VzVzg{r,-n.,uj), (101) 

where g{r,n,uj) is the scalar Green function of the Helmholtz equation (j98p . satis- 
fying the radiation condition at infinity and the finiteness condition on the axis of 
the cylinder. Such a function is known to be given by the following series expansion 

X / dhIp{vr)Kp{vRcN) e"'^'-'"\ r < Rcn, 
Jc 

where Ip and Kp are the modified cylindric Bessel functions, v = v{h, uj) = Vh"^ — k"^, 
and we used the property (j82l) to go from the discrete variable n to the corresponding 
continuous variable. The integration contour C goes along the real axis of the 
complex plane and envelopes the branch points it/c of the integrand from below and 
from above, respectively. For r > Rcn, the function g(r, n, w) is obtained from 
Eq. (I102p by means of a simple symbol replacement Ip o Kp in the integrand. 

The scalar function (jl02p is to be imposed the boundary conditions on the CN 
surface. To derive them, we represent the classical electric and magnetic field com- 
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ponents in terms of the EM field Green tensor as follows 

E^{r,uj) = ik^Gaz{r,n,u), (103) 

H^{r,io) = -^ ^ e^^^V^E^{r,uj). (104) 

These are valid for r 7^ n under the Coulomb-gauge condition. The boundary condi- 
tions are then obtained from the standard requirements that the tangential electric 
field components be continuous across the surface, and the tangential magnetic 
field components be discontinuous by an amount proportional to the free surface 
current density, which we approximate here by the (strongest) axial component, 
<^zz{RcN ,^)^ of the nanotube's surface conductivity. Under this approximation, 
one has 

E,\+-E,\.=E^\+-E^\- = Q, (105) 

FJ+-^J_=0, (106) 
47r 

H^\+ - H^\. = — azz{io)EjR^^, (107) 

where it stand for r = Rqn^^ with the positive infinitesimal e. In view of Eqs. ()103p . 
P04|) and pOO|) . the boundary conditions above result in the following two boundary 
conditions for the function (jl02p 



5l+-5l_=0, (108) 



dg 

dr 



dg 

dr 



^2 



+ ^ MUc^- (109) 



We see that Eq. (llOSp is satisfied identically. Eq. (I109P yields the Wronskian of 
modified Bessel functions on the left, W[Lp{x),K,p{x)\=Ip{x)K'p{x) — Kp[x)I'p{x) = 
— 1/x, which brings us to the equation 

1 _4^^^-('^)^2^^(^^^^)^^(^^^^). (110) 



RcN w 

This is nothing but the dispersion relation which determines the radial wave num- 
bers, /i, of the CN surface EM modes with given p and u. Since we are interested 
here in the EM field Green tensor on the CN surface [see Eq. (j93p ]. not in particular 
surface EM modes, we substitute Ip{vRcN)Kp{vRcN) from Eq. (jllOp into Eq. (jl02p 
with r = RcN- This allows us to obtain the scalar Green function of interest with 
the boundary conditions (jlOSp and (jlOOp taken into account. We have 

c/(R,n,a;) = -— dh—, (111) 



8iT'^azz{i^)RcN Jc k'^ -h? ' 
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where R = {Rcn,^,z} is an arbitrary point of the cyhndrical surface. Using 
further the residue theorem to calculate the contour integral, we arrive at the final 
expression of the form 



which yields 



,(R,n,.) = - ^y^y,7"^ e-l^-l/^ (112) 



^C{n,n,uj)=0, (113) 
llG,,(R,n,t^) =<7(R,n,w), (114) 

in view of Eqs. ifTOO]) and (fTOTD . 

The fact that the transverse Green function ()113p identically equals zero on the 
CN surface is related to the absence of the skin layer in the model of the infinitely 
thin cylinder (see, e.g., Ref. [HSj). In this model, the transverse Green function 
is only non-zero in the near-surface area where the exciton wave function goes to 
zero. Thus, only longitudinally polarized EM modes with the Green function ()114p 
contribute to the exciton surface EM field interaction on the nanotube surface. 



Appendix C 

Diagonalization of the Hamiltonian (fll)— (1131) 

We start with the transformation of the total Hamiltonian (fT]l~ (113p to the k-repre- 
sentation using Eqs. (I5|) and ([9]). The unperturbed part presents no difficulties. Spe- 
cial care should be given to the interaction matrix element gj^^(n, m, w) in Eq. (I13p . 
In view of Eqs. (|114p . (|112p and (l3|), one has explicitly 

gf )(k,k',a;) = ^ J^gf )(n,m,c^)e-*''-"+*^'- (115) 



n,m 



lUJ 



V^7r^Re(j^^(w) di i?- 



— V Jo 



27rca,,(a;)i?ciV N " NSl 

Jo J ~oo 

where we have also taken into account the fact that the dipole matrix element 
(d;^)^ = (OKdn)^!/) is the same for all the lattice sites on the CN surface in view of 
their equivalence. As a consequence, (d;^)^ = dl/N with dl = X^n(0|(dn)^|/). 
The integral over ip in Eq. (jll5p is taken in a standard way to yield 

d^nd^m5{^n- ^m) e-^'=^'^"+^^>™ = 27r5fc^fc^. (116) 
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The integration over z is performed by first writing the integral in the form 



Um 



L/2 j-L/2 
dZn I dz. 



'-L/2 J-L/2 

(L being the CN length), then dividing it into two parts by means of the equation 



Ju)\z„-Zm.\/c 



and finally by taking simple exponential integrals with allowance made for the for- 
mula 

2sin[L(/c^- A;;)/2] 



lim 



L->oo L{kz — k'^) 
After some simple algebra we obtain the result 



dzndzme''^^'--''^^/^-'''^^"-^''''-^"^ = lim L^\l- — 

L— >oo I L 



2iu/c 



(117) 



In view of Eqs. (jll6p and (|117p . the function (|115p takes the form 



gf^(k,k',..) 



di \J TrSohuj Re 



lim si — 

{2iT)'^cazz{'^)RcN L^ooi L[k1—{uj/cY^ 



2iuj/c 



Skk'. (118) 



We have taken into account here that 5^^^' ^k^k 

\2 _ 1 U'^^\2 



(5kk/, as well as the fact that 
{RcnL/NSq)'^ = l/(27r)^. This can be further simplified by noticing that only 
absolute value squared of the interaction matrix element matters in calculations of 
observables. We then have 



2iuj/c 



L\ 



a 



1 + ^-1 + 



a 



V? + cP' 



with u = {ckz/ijj)'^— 1, and a = {2c/ Luj)"^ being the small parameter which tends to 
zero as L — )• oo. Using further the formula (see, e.g., Ref. [73] ) 

5{u) = — lim 



and the basic properties of the 5-function, we arrive at 



lim 



1 



2iuj/c 



L 



[u/cf 



1 + 



[6{u} + ckz) + S{uj - ckz)] 



We also have 



Y^Re azz{' 



OJ 



Re 



crzz[U} 



Cfzzil^) ' 



(119) 



(120) 



46 



I.V.Bondarev, L.M.Woods, and A.Popescu 



Equation (jllSp . in view of Eqs. (|119p and p20p . is rewritten effectively as follows 

gf)(k,k',a;) = ±zI?/H4k' (121) 

with 

Df{uj) = (27r)^cRcN V ^ ~ ^ ' ^ ^ 

In terms of the simplified interaction matrix element (|12ip , the k-representation 
of the Hamiltonian (fTI)- (|13p takes the following (symmetrized) form 



(123) 



where 



(124) 



+ J^du;nLO [/t(k, a;)/(k, w) + /t(-k, u;)f{-k, u) 
+ ^^/H (4,/ + ^-kj) [/(k,^) - fH-K 



CO 



+ /i.e. 



with Df{(jj) given by Eq. (jl22p . To diagonalize this Hamiltonian, we follow Bo- 
goliubov's canonical transformation technique (see, e.g., Ref. [7^). The canonical 
transformation of the exciton and photon operators is of the form 



i?k,/ = [^^(k,a;/)e;(k) +t;^(k,a;/)et(_k)" 
/(k,..) = E [n;(k,a;)e;(k) + ^;;(k,L.)|t(-k) 



(125) 
(126) 



^l=l,2 



where the new operators, |^(k) and |Jj(k) = [|^(k)]"l', annihilate and create, respec- 
tively, the coupled exciton-photon excitations of branch fj, on the nanotube surface. 
They satisfy the bosonic commutation relations of the form 



(127) 



which, along with the reversibility requirement of Eqs. ()125p and (jl26p . impose the 
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following constraints on the transformation functions and 
K(k,w/)n^'(k,a;/) - ?;^(k,a;/)t;*,(k, w/)] 

/ 

/•oo 

+ duj [u^{k,u})u*,{k,u}) - v* {k,uj)v^/{k,uj)] = (5^^/, 
Jo 



Here, the first equation guarantees the fulfilment of the commutation relations ()127p . 
whereas the second and the third ensure that Eqs. (|125p and (|126p are inverted to 
yield ^/^(k) as given by Eq. (|15p . Other possible combinations of the transformation 
functions are identically equal to zero. 

The proper transformation functions that diagonalize the Hamiltonian ()124p to 
bring it to the form (|14p . are determined by the identity 



(128) 



Putting Eqs. (|15p and (|124p into Eq. (jl28p and using the bosonic commutation 
relations for the exciton and photon operators on the right, one obtains (k-argument 
is omitted for brevity) 



{hw^ - Ef) u*{ujf) = -i dujDf{uj) [u^{uj) - v*{uj)\ , 

Jo 

/•oo 

{huj^ + Ef) Vf,{u}f) =i duj Df{uj) [u^{uj) - v*{uj)\ , 
Jo 

h (w^ - uj) u^{uj) = iY Dfiuj) [ul^iujf) + v^iuf)] , 
/ 

h {u}f, + uj) V* (cj) = Df{u;) [u*^{ujf) + Vf,{ujf)] . 
f 

These simultaneous equations define the complex transformation functions and 
uniquely. They also define the dispersion relation (the energies hw^, fj, = 1,2) 
of the coupled exciton-photon (or exciton-plasmon, to be exact) excitations on the 
nanotube surface. Substituting and from the third and forth equations into 
the first one, one has 



fiu},, - Ef 



AE 



f 



du! 



uj\Df {uj)\ 



2 1 
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whereby, since the functions u* are non-zero, the dispersion relation we are interested 
in becomes 

-E)- AEf Hdu 1^^^ = . (129) 

Jo K^i - ^ ) 

The energy Eq of the ground state of the coupled exciton-plasmon excitations is 
found by plugging Eq. (fTSj) into Eq. (fT4|) and comparing the result with Eqs. (1123P 
and (fT2iD . This yields 



EQ = -Y,^^.{^) 

k, n=l,2 



f 



POO 

\Vf_,(k,ujf)\'^ + / duj\v^{k,uj 
Jo 



Using further Df{uj) as explicitly given by Eq. (jl22p . the dispersion relation 
(I129P is rewritten as follows 

Here we have taken into account the general property (Tzz{^) = w), which 

originates from the time-reversal symmetry requirement, in the second term on the 
right hand side. This term comes from the two delta functions in |Z)j(a;)p, and 
describes the contribution of the spatial dispersion (wave- vector dependence) to the 
formation of the exciton-plasmons. We neglect this term in what follows because 
the spatial dispersion is neglected in the nanotube's axial surface conductivity in our 
model, and, secondly, because it is seen to be very small for not too large excitonic 
wave vectors. Thus, converting to the dimensionless variables (I17p . we arrive at the 
dispersion relation (jl6p with the exciton spontaneous decay (recombination) rate 
and the plasmon DOS given by Eqs. (fT8]l and (fT9]) . respectively. 

Lastly, bearing in mind that the delta functions in \Df{ijj)\'^ are responsible 
for the spatial dispersion which we neglect in our model, and therefore dropping 
them out from the squared interaction matrix element (jl2ip . we arrive at the prop- 
erty (US]). 



Appendix D 

Effective longitudinal potential in the presence 
of the perpendicular electrostatic field 

Here we analyze the set of equations (|27p - (j29p . and show that the attractive cusp- 
type cutoff potential ((32]) with the field dependent cutoff parameter ((33]) is a uni- 
formly valid approximation for the effective electron-hole Coulomb interaction po- 
tential (130]) in the exciton binding energy equation (|29]) . 
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We rewrite Eqs. (j271) and (j28|) in the form of a single equation as follows 



dip^ 



+ +pcosip]^{ip) = 0. (130) 



Here, if = ipe,h, ip = i'e.h, q = RcNs/^m^J^e^/h, and p = ±2eme^hR^cj^F /h^ with 
the (+)-sign to be taken for the electron and the (-)-sign to be taken for the hole. We 
are interested in the solutions to Eq. (jl30p which satisfy the 27r-periodicity condition 
tp[ip) = ip({p + 2tt). The change of variable ip = 2t transfers this equation to the 
well known Mathieu's equation (see, e.g., Refs. [1191 1114] ). reducing the solution's 
period by the factor of two. The exact solutions of interest are, therefore, given 
by the odd Mathieu functions se2m+2{t = Lp/2) with the eigen values ?)2m+2i where 
m is a nonnegative integer (notations of Ref. |119j ). These are the solutions to 
the Sturm-Liouville problem with boundary conditions on functions, not on their 
derivatives. 

It is easier to estimate the z-dependence of the potential (j30p if the functions 
'4'e,hiVe,h) are known explicitly. So, we do solve Eq. (|13U|) using the second order 
perturbation theory in the external field (the term p cos (p) . The second order field 
corrections are also of practical importance in the most of experimental applications. 

The unperturbed problem yields the two linearly independent normalized eigen 
functions and the eigen values as follows 

(0) exp(ifcijV) . RcN fZ (oy 

with j being a nonnegative integer. The energies are doubly degenerate 

with the exception of ££^^^(0) = 0, which we will discard since it results in the 
zero unperturbed band gap according to Eq. ([8]). The perturbation pcosip does 
not lift the degeneracy of the unperturbed states. Therefore, we use the standard 
nondegenerate perturbation theory with the basis wave functions set above (plus 
sign selected for definiteness) to calculate the energies and the wave functions to the 
second order in perturbation. The standard procedure (see, e.g., Ref. [95j) yields 

'^^^^'^^^^'^^ = [ [U-ir-Pr + [(, + l)2_,2]2 j 2h^ P ^,e,Hi^e,H) 



' ^(j-2)<\,,,(v.eA) i^flleA^e,k) 



^(j-2)^(j-3)V>fJ,,^,(^e,,) i^%e,h{Ve,H) \ K,H^^R'cN ^2 
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Here, j is a positive integer, and the theta-functions ensure that j = 1 is the ground 
state of the system. The corresponding energies are as foUows 



with Wj given by Eq. (|Mj) . thus, according to Eq. ([8]), resulting in the nanotube's 
band gap as given by Eq. ([57|) . 

From Eq. ()132p , in view of Eq. (jlSip , we have the foUowing to the second order 
in the field 



, \ ^ F^f^ AT IV j , . 

1 - 2(m/,cos(^/, - meCOS99e) — ^t^-^ (134) 



+2 (m^ cos 2(/5/, + cos 2ipe) " ^M^ex cos (^e COS (^/j J^r"^^ 

where 

^?a-2) r ??(i-3) , 1 



Plugging Eqs. (jl34p and ()3ip into Eq. pop and noticing that the integrals involv- 
ing linear combinations of the cosine-functions are strongly suppressed due to the 
integration over the cosine period, and are therefore negligible compared to the one 
involving the quadratic cosine-combination, we obtain 

1/(1 l-2cosy.cosytAj(F) 



with Aj(F) given by Eq. 

The next step is to perform the double integration in Eq. (jl35p . We have to 
evaluate the two double integrals. They are 

and 

"^J f ^"" d^h cos V?e COS iph , 

/o ^7o{z2 + 4ii2^^sin2[(^,-(^,)/2]}V2 • ^ 
We first notice that both Ii and I2 can be equivalently rewritten as follows 

/'27r /'27r p27r r^e. 

/ dLpei diph--- = 2 / d(/9e/ rfv'/i--- (138) 
Jo Jo Jo Jo 
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due to the symmetry of the integrands with respect to the (93e = 93/i)-hne. Using this 
property, we substitute iph with the new variable t = sm[{ipe — ^h)/"^] in Eqs. ()136p 
and (|137p . This, after simphfications, yields 



and 



■4.1 d^, 

Jo 

2n 



sin((/Pe/2) 



dt 



/2 = 4 / dife cos (ft 



[(l-t2)(,2+4^2^^i2)]l/2 
sin(^e/2) rfi(l_2t2) 



(139) 



(140) 



/O [(l-t2)(z2 + 4^2^^i2)]l/2 • 

Here, the inner integrals are reduced to the incomplete elliptical integrals of the first 
and second kinds (see, e.g., Ref. |114j ). 

We continue the evaluation of Eqs. (|139p and ()140p by expanding the denom- 
inators of the integrands in series at large and small \z\ as compared to the CN 
diameter 2Rcn- One has 



1 



1 



(z2+4i?^^t2)l/2 |, 

for \z\/2RcN > 1, and 
Jo 



1 _ i ( ^^CNt y 3 f 2RcNt y 5_ (2RcNt\ 
2\ \z\ J 8 V kl / 16 V \A J 



+ ... 



lim 



sm((/Pe/2) j^^^ 



[(1 _ t2)(^2 + ARl^t^)Y/^ 2RcN (kl/2Kciv)^0 7|.|/2RciV tVl^ 

for \z\/2RcN ^ 1 [f{t) is a polynomial function]. Using these in Eqs. ()139p and 
(jl40p . we arrive at 



h 



RcN 



and 



47r' 



RcN 



< 1 



> 1 



1 + 



]^ _ 3 / 2RcN 



2" 



27^ 



' 2i?c7V 
> 1 



< 1 



Plugging these Ii and /2 into Eq. (jl35p and retaining only leading expansion terms 
yields 



f e^[l-A,(F)] i^^l^^i 



^eff(^) 



neRcN 



e\z\ ' 272(7 AT 



> 1 



|z| y 2i2c7V 



< 1 



(141) 
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Figure 14: The dimensionless function (|142p with the zero-field cutoff parameter 
(fTi3]) . See text for details. 

We see from Eq. (jl4ip that, to the leading order in the series expansion parame- 
ter, the perpendicular electrostatic field does not affect the longitudinal electron-hole 
Coulomb potential at large distances \z\ » 2Rcn, as one would expect. At short 
distances \z\ <C 2Rcn the situation is different, however. The potential decreases 
logarithmically with the field dependent amplitude as \z\ goes down. The amplitude 
of the potential decreases quadratically as the field increases [see Eq. (I34p ]. thereby 
slowing down the potential fall-off with decreasing \z\, or, in other words, making the 
potential shallower as the field increases. Such a behavior can be uniformly approx- 
imated for all \z\ by an appropriately chosen attractive cusp-type cutoff potential 
with the field dependent cutoff parameter. Indeed, consider the dimensionless func- 
tion f{y) = — 2i?c'Are V^fj/e^ of the dimensionless variable y = \z\/2Rcn- Then, 
according to Eq. p4ip . one has 



fiy) = 

Now introduce the function 



<I>i(y) = |[l-A,(F)]ln(|), 0<y«l 



My) = ^, y > 1 



Hy) = ^— (142) 

y + yo 

with the cutoff parameter yo selected in such a way as to satisfy the condition 
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$(1) = [$i(l) + «>2(l)]/2. This yields 

^° - vr + 21n2[l-A,(F)] ' ^^^^^ 

Figure [T^ shows the zero-field behavior of the $(y) function as compared to the 
corresponding and ^*2(y) functions. We see that ^{y) gradually approaches 

^2iy) = 1/y for increasing y > 1. For decreasing y < 1, on the other hand, 
is very close to the logarithmic behavior as given by <&i(y), with the exception 
that there is no divergence at y ~ due to the presence of the cutoff. The cutoff 
parameter (jl43p is field dependent, decreasing as the field grows, which is consistent 
with the behavior of the original potential (jl4ip . Multiplying Eq. (I142p by the 
dimensional factor —e'^/2RcN^ and putting y = \z\/2Rcn-, we obtain the attractive 
longitudinal cusp-type cutoff potential ()32p we build our analysis on in this paper. 
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